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SECTION  I 


INTRODUCTION 

Due  to  in-flight  emergencies  or  other  unforseen  events  it  occa¬ 
sionally  becomes  necessary  for  Air  Force  aircraft  to  discharge 
unburned  fuel  directly  into  the  atmosphere  while  airborne.  The 
primary  reason  for  Jettisoning  fuel  is  to  reduce  the  aircraft’s 
gross  weight  to  facilitate  a  safe  landing.  As  a  result  of  a 
growing  awareness,  of  fuel  Jettisoning ' s  potential  for  environmen 
tai  degradation,  the  Air  Force  initiated  a  study  in  1972  to  deter 
mine  the  nature,  the  extent  and  the  environmental  affects  of  fuel 
Jettisoning  by  Air  Force  aircraft.  The  overall  results  of  that 
study  are  reported  in  Reference  1.  The  present  report  documents 
the  experimental  and  modelling  efforts  which  were  undertaken  to 
determine  the  physical  fate  of  the  Jettisoned  fuel  in  support  of 
one  aspect  of  that  study. 

Jet  feel  ,  when  jettisoned  from  an  aircraft,  readily  breaks  up 
into  small  droplets  and  begins  to  evaporate.  The  fuel  vapor  and 
droplets  are  subject  to  entrainment  in  the  aircraft  wake,  disper¬ 
sion  by  atmospheric  turbulence,  and  (in  the  case  of  droplets) 
gravitational  settling.  The  principal  question  to  be  answered  by 
the  work  described  in  this  report  was  what  fraction  of  the 
jettisoned  fuel  can  be  expected  to  reach  the  ground  in  liquid 
form.  The  environmental  affects  due  to  liquid  fuel  contamination 
of  the  ground  are  very  different  from  those  associated  with  the 
airborne  fuel  vapors.  Therefore,  a  knowledge  of  the  fraction  of 
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fuel  In  each  form  Is  essential  for  estimation  of  the  .erall 
impact  of  fuel  Jettisoning. 


Section  II  describes  a  computer  model  which  was  developed  to  simuj 
late  the  evaporation  and  free-fall  of  fuel  droplets  in  the 
atmosphere.  In  order  to  apply  this  model  to  fuel  jettisoning  it 
"is  necessary  to  know  the  distribution  of  droplet  sizes  which  is 
^produced  by  the  jettisoning  process;  the  experimental  study  which 
was  performed  to  determine  this  information  is  described  in 
Section  III.  As  noted  there  it  was  necessary  to  use  the  droplet 
evaporation  model  to  interpret  the  experimental  findings.  During 
this  experiment,  sampling  was  also  performed  at  ground  level  to 


determine  whether  the  jettisoned  fuel  reached  the  ground  In  signi¬ 


ficant  concentrations.  The  results  of  the  ground-sampling  are 
also  presented  in  this  section.  '  “ 


In  Section  IV  the  predictions  of  the  droplet  model  concerning  the 
evaporation  and  groundfall  of  jettisoned  fuel  are  presented  along 
with  an  estimate  of  the  atmospheric  dispersion.  A  detailed 
description  of  the  droplet  model,  together  with  a  program  listing 
and  sample  output,  is  provided  in  the  appendices. 


SECTION  II 


FUEL  DROPLET  FREE-FALL  AND  EVAPORATION  MODELLING 

The  evaporation  of  a  fuel  droplet  is  complicated  by  the  fact  that 
as  the  droplet  evaporates  its  composition  changes.  The  lighter, 

-  S' 

more  volatile  components  are  stripped  away  preferentially  while 
the  denser,  slower-evaporating  components  remain.  Also,  as  the 
droplet  evaporates  it  becomes  smaller,  reducing  the  surface  area 

.  available  for  evaporation.  Therefore  the  evaporation  rate 

decreases  markedly  at  first,  but  later  levels  off.  A  further 
complication  arises  when  the  droplet  is  allowed  to  fall.  As  it 
falls  the  droplet  experiences  changing  atmospheric  temperature, 
pressure,  and  viscosity,  since  these  properties  are  a  function  of 
altitude.  These  changing  atmospheric  properties  modify  the 
droplet’s  evaporation  rate  as  well  as  its  rate  of  fall.  In 
addition,  the  rate  of  fall  varies  with  the  droplet’s  size  and 
density,  which  in  turn  depend  on  its  evaporative  history.  As  a 
result,  the  evaporation  and  free-fall  of  the  droplet  are  inter¬ 
dependent  processes. 

2.1  DESCRIPTION  OF  THE  MODEL 

The  model  which  was  developed  to  simulate  the  simultaneous  evapor 
ation  and  free-fall  of  a  fuel  droplet  is  based  on  the  pioneering 
work  of  Lowell  (References  2  and  3).  Lowell's  method  has  been 
refined  and  extended  chiefly  by  the  incorporation  of  experimental 
data  provided  by  the  Arnold  Engineering  Development  Center  (AEDC) 
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and  by  the  use  of  a  more  detailed  fuel  composition.  The  data 
obtained  by  AEDC  (Reference  4)  includes  measurement  of  the  drag 
coefficients  and  evaporation  rates  of  JP—  4  fuel  droplets  as  well 
'as  the  physical  properties  of  JP-4  as  a  function  of  temperature. 
Data  on  the  composition  of  JP-4  provided  by  the  Air  Force  Aero- 
Propulsion  Laboratory  was  used  to  prepare  a  theoretical  33  com¬ 
ponent  ’’synthetic  JP-4’’  mixture  as  shown  in  Table  1.  The 
calculated  physical  properties  of  this  mixture  (density,  average 
carbon  number,  and  Reid  vapor  pressure)  egree  well  with  typical 
analyses  of  JP-4.  Each  component  in  the  synthetic  mixture  repre¬ 
sents  a  class  of  compounds  in  real  JP-4.  Lowell  used  ten  equal 
volume  distillate  fractions  to  represent  JP-4;  the  more  detailed 
composition  allows  modelling  of  the  droplet  to  continue  until  vir¬ 
tually  all  ( 99 • 9% )  of  the  mass  is  gone.  Additional  refinements 
include  a  complete  energy  balance,  and  an  initial  droplet  tem¬ 
perature  estimate.  The  effect  of  these  two  calculations  is 
assessed  in  Appendix  A. 

Essentially,  the  model  breaks  up  a  droplet’s  fall  into  a  series  of 
small  time  intervals.  The  distance  of  fall  during  each  interval 
is  calculated  assuming  the  droplet  is  falling  at  the  terminal 
velocity  for  its  current  diameter,  density  and  altitude.  Loss  of 
mass  through  evaporation  is  calculated  assuming  Raoult’s  law; 
that  is,  each  component  evaporates  independently.  An  energy  bal¬ 
ance  routine  adjusts  the  droplet  temperature  to  allow  for  evapora¬ 
tive  cooling,  radiation,  conduction,  and  insolation  effects.  The 
new  droplet  composition,  mass  and  altitude  are  used  as  initial 
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TABLE  1.  SYNTHETIC  JP-4 


C ompound 


Volume 

Percent 


Molecular 

Weight 


Bolling 
Point (C ' 


Density 

(g/ml) 


iso-pentane 

iso-hexane 

cyclohexane 

benzene 

\  3-methylhexane 

methylcyclohexane 

toluene 

•  . 4-methylheptane 

.  cis-1 , 4-diethylcyclohexane 
m-xylene 
4-methyloctane 
lsop ropy 1 eye lohe  xane 

1- ethyl-2-methylbenzene 
2, 7-dimethyloctane 
p-menthane  (cis) 
p-cymene 

:  napthal ene 
undecane 

3-methylbutylcyclohexane 

2- methylenedecalin  (trans) 
_ l-butyl-3-methylbenzene 

1 -methylnapthalene 
dodecane 

3- e thy lbutyl cyclohexane 

1 . 3. 5 - triethylbenzene 
2, 3-dimethylnapthalene 

— tridecane 

3-isoprop.y  L  butyl  eye  lohe  xane 

3.5- di.  thyl-l-propylbenzene 
tetradecane 
pentadecane 
perhydrophenanthrene 
pyrene 


;  3.9 

72.2 

27.9 

.62 

8.1 

86.2 

60.2 

.66 

2.1 

84.2 

80.7 

.78 

0.3 

78.1 

80.1 

.88 

9.4 

100.2 

91.8 

.69 

7.1 

9  8.2 

100.9 

.77 

••  0.7 

r  92.1 

ml  10.8 

-.87 

10.1 

114.2 

117.7 

--.70 

7.4 

112.2 

"  124.3 

.78 

1.6 

106.2 

139.1 

.87 

9.1 

128.3 

142.4 

.  .72 

4.3 

126.2 

154.5 

.80 

2.4 

120.2 

165.2 

.88 

7.3 

142.3 

159.6 

.72 

3.7 

140.3 

170.9 

.80 

1.8 

134.2 

177.1 

.86 

0.2 

128.2 

-217.9 

1.03 

4.8 

156.3 

195*9 

.74 

2.5 

154.3 

196.5 

.30 

3.4 

150.3 

201.0 

.89 

—1.1 

148.2 

205.0 

r  .86 

0.2 

142.2 

244.6 

1.02 

2.8 

170.3 

216.3 

•  75 

1.2 

168.3 

211.0* 

.80* 

0.5 

162.3 

216.0 

.86 

0.2 

156.2 

268.0 

1.00 

1.1 

184.4 

235.4 

•  7  6 

0.4 

.  182.4 

225.0* 

.80* 

0.1 

176.3 

234.0* 

.87* 

0.2 

198.4 

253.7 

.75 

0.1 

212.4 

270.6 

.77 

1.8 

192.4 

290.0* 

.94 

0.1 

202.3 

393.0 

1.27 

Density  =  0.75  g/ml  NOTE:  Olefins  can  be  con 

Average  Carbon  #  =  8.7  sidered  present  In  the 

Reid  vapor  pressure  (100.F)  =2.5  psla  paraffinic  fractions  at  a 
Aromatics  =  9 • 2%  level  of  1$. 

Napthalenes  =  0.7 % 

^estimated 
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conditions  for  the  next  interval.  This  stepwise  approximation 
continues  until  the  droplet  impacts  on  the  ground  or  loses  99*9 
percent  of  its  initial  mass.  The  initial  conditions  which  must  be 
known  are  the  droplet’s  original  composition,  altitude  and  diam¬ 
eter,  the  temperature  at  local  ground  level,  and  the  aircraft's 
air  speed.  The  initial  droplet  temperature  is  then  taken  as  the 

-corresponding  stagnation  temperature,  assuming  equilibration  of _ 

the  fuel  tanks  with  the  skin  of  the  aircraft.  In  the  early  inter¬ 
vals  the  droplet  is  then  allow  to  cool  through  evaporation  until 


an  energy  balance  is  achieved. 

A  more  detailed  description  of  the  model  is  presented  in  Appendix 

B,  and  a  listing  of  the  computor  program  is  included  as  Appendix 

C.  The  program  reports  the  droplet's  status  at  the  end  of  each 
time  interval,  or  a  selected  period  of  time.  Output  for  a  single¬ 
droplet  case  includes  the  elapsed  time  since  release;  the  droplet 
altitude,  velocity,  diameter  and  fraction  of  initial  mass 
remaining;  and  the  fractional  mass  of  each  component.  Sample 
output  for  several  single-droplet  cases  as  well  as  for  a  distri¬ 
bution  of  droplets  is  provided  in  Appendix  D. 

2.2  SINGLE-DROPLET  RESULTS 

Figure  1  shows  the  predicted  evaporation  and  free-fall  of  a  270- 
mlcron  droplet  released  at  1500  meters  for  three  different  ground- 
level  temperatures.  Initially  evaporation  dominates,  and  in  the 
first  few  minutes  the  droplet  loses  from  60  to  90  percent  of  its 
mass,  depending  on  the  temperature.  Thereafter  evaporation  be¬ 
comes  less  important,  and  the  droplet  falls  without  a  substantial 
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change  in  mas3.  Ilhe  duration  of  the  evaporation  phase  for  all 
three  temperatures  is  around  five  minutes,  while  the  free-fall 
phase  ranges  from  l;50  minutes  at  -20°C  to  over  11  hours  at  +20°C. 

At  the  higher  temperature  the  droplet  evaporates  more  completely 
in  the  first  phase;  and  therefore,  being  smaller,  it  falls  more 
slowly  during  the  second  phase.  The  dots  on  the  curve  for  0°C  show 
the  actual  intervals  used  by  the  model  in  performing  its  stepwise 
simulation.  The  model  maintains  roughly  constant  increments  of 
mass  and  then  altitude  by  increasing  the  duration  of  the  intervals 
from  less  than  a  tenth  of  a  second  initially  to  more  than  30  rnin- 
iites  at  the  end. 

The  effect  of  a  droplet's  initial  size  is  shown  in  Figure  2  for 
an  initial  altitude  of  1500  meters  and  a  ground-level  temperature 
of  G°C.  On  a  relative  basis,  a  100-micron  droplet  appears  to  evap¬ 
orate  more  quickly  than  a  500-mlcron  droplet.  However,  this  is 
only  an  artifact,  caused  by  the  fact  that  the  mass  of  a  droplet  is 
proportional  to  its  diameter  cubed  while  the  evaporation  rate  is 
proportional  to  the  diameter  squared.  On  an  absolute  basis,  the 
larger  droplet  is  actually  evaporating  much  faster  than  the 
smaller  droplet,  as  can  be  seen  from  the  more  rapid  change  in 
diameter.  Due  to  its  larger  size,  the  500-micron  droplet  also 
falls  much  faster  than  the  100-micron  droplet. 

As  mentioned  earlier,  the  more  volatile  components  in  the  fuel 
evaporate  preferentially.  Figure  3  demonstrates  this  effect  for 
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a  270-micron  droplet  released  at  1500  meters  when  the  ground-level 
temperature  is  -20°C.  The  lighter  compounds  (containing  less  than 
eight  carbon  atoms)  are  lost  almost  immediately,  while  the  heavier 
components  (containing  more  than  twelve  carbons)  hardly  evaporate 
at  all.  This  results  in  the  different  fuel  components  being 
distributed  vertically:  the  lighter  compounds  are  left  at  the 
initial  release  height,  the  heavier  compounds  reach  the  ground, 
and  the  intermediate  compounds  are  spread  out  in  between. 

Finally,  Figure  4  examines  the  related  question  of  whether  the 
evaporation  of  a  complex  mixture  can  be  approximated  by  a  single 
compound.  In  Figure  3  the  curves  C6  through  Cl4  represented  dif¬ 
ferent  components  in  the  same  droplet.  In  Figure  4,  however, 
curves  C6  through  Cl  4  represent  different  droplets,  each  of  which 
contain  only  one  component.  These  curves  show  that  no  single- 
component  droplet  matches  the  behavior  of  a  droplet  of  JP-4. 

During  the  early  stages  of  evaporation,  a  JP-4  droplet  resembles 
hexane  or  octane.  But  then  the  JP-4  droplet  evaporation  slows 
down  until  it  eventually  parallels  that  of  a  dodecane  or  tetrade- 
cane  droplet.  From  inspection  of  the  results  for  0°C  and  -20°C  it 
can  be  seen  that  no  single-component  droplet  could  be  used  to 
predict  the  fraction  of  JP-4  reaching  the  ground  under  different 
conditions.  For  the  same  reason,  the  ability  of  a  mixture  of 
several  compounds  to  represent  JP-4  is  limited  when  the  number  of 
compounds  in  the  mixture  is  small.  For  example,  a  mixture  of  ten 
equal  fractions  could  not  simulate  a  JP-4  droplet  beyond  90  per¬ 
cent  evaporation,  because  at  that  point  the  simulated  droplet 
would  essentially  consist  of  just  one  component. 


2.3  COMPARISON  WITH  EXPERIMENTAL  DATA 


One  of  the  experiments  carried  out  by  AEDC  in  support  of  this  fuel 
jettisoning  study  consisted  of  suspending  small  droplets  of  JP-4 
from  a  glass  filament  and  measuring  their  rates  of  evaporation 
when  subjected  to  a  controlled  airflow  (Reference  4).  The  results 
obtained  for  several  droplets  are  shown  as  points  in  Figure  5. 

-These  data  were  obtained  at  a  temperature  of  20°C  and  with  an 
airflow  of  3  meters  per  second,  which  is  close  to  the  terminal 
velocity  of  the  size  droplets  studied.  The  curves  in  Figure  5 
show  the  predictions  of  the  fuel  droplet  model  under  the  same 
conditions.  The  agreement  is  excellent. 

The  AEDC  effort  also  attempted  to  measure  evaporation  rates  in 
cooled  airflows  (-10°  to  -20°C).  Unfortunately,  water  conden¬ 
sation  at  these  temperatures  led  to  severe  experimental  problems, 
bringing  the  validity  of  the  few  results  obtained  into  question. 
The  reported  evaporation  rates  are  extremely  low,  roughly  a  factor 
of  five  lower  than  the  fuel  droplet  model  would  predict.  Since 
the  temperature  dependent  parameters  in  the  fuel  droplet  model 
have  been  checked  against  experimental  values  (see  Appendix  B), 
the  disagreement  solely  with  AEDC's  low-temperature  data  indicates 
that  the  low-temperature  experiments  may  have  been  compromised  by 
the  water  condensation  problems.  However,  these  problems  were  not 
encountered  during  the  tests  at  +20°C,  so  the  use  of  the  room- 
temperature  results  to  validate  the  fuel  droplet  model  appears  to 
be  justified. 


SECTION  III 


EXPERIMENTAL  FUEL  JETTISONING  STUDY 
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The  fuel  droplet  model  described  in  the  previous  section  begins 
with  an  initial  droplet  diameter,  altitude  and  temperature  and 
predicts  the  ensuing  history  of  the  droplet.  Therefore  the 
application  of  this  model  to  the  question  of  fuel  jettisoning 
■  requires  a  knowledge  of  the  droplet  sizes  Involved.  It  can  be 

assumed  that  in  the  process  of  being  discharged  through  an  orifice 
into  a  high  velocity  airstream,  the  fuel  would  be  atomized  into  a 
...distribution  of  droplet  sizes.  However,  the  droplets  produced 
could  range  from  only  a  few  microns  in  diameter  (the  size  that 
scatter  light)  up  to  several  thousand  microns  in  diameter  (the 
size  of  raindrops).  The  actual  droplet  size  distribution  produced 
by  the  jettisoning  process  cannot  be  predicted  from  theory. 
Therefore  it  must  be  determined  experimentally. 

3.1  BACKGROUND 

An  initial  attempt  to  determine  the  droplet  size  distribution 
produced  by  fuel  jettisoning  was  performed  in  1973  (Reference  5). 
This  study,  which  was  conducted  at  the  AEDC,  simulated  the  fuel 
jettisoning  process  by  injecting  JP-4  into  a  high-velocity 
airstream  in  a  wind  tunnel.  Airflows  of  100  to  200  meters  per 
second  were  used,  representing  typical  aircraft  airspeeds  during 
jettisoning.  The  droplets  produced  were  then  measured  using 
holography.  The  observed  droplet  size  distribution  is  shown  in 
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Figure  6  ( "AEDC-experimental" ) .  All  of  the  droplets  detected  were 
.  less  than  100  microns  in  diameter,  with  a  mass  median  diameter  of 

only  32  microns.  However,  due  to  the  constraints  imposed  by  the. _ 

size  of  the  wind  tunnel  these  tests  were  conducted  with  fuel  flow 
rates  of  only  0.02  to  2.2  kilograms  per  second,  compared  to  typi¬ 
cal  aircraft  jettisoning  rates  of  5  to  50  kilograms  per  second, 
j  .Therefore,  it  was  not  certain  that  these  results  were  represen¬ 

tative  of  full-scale  fuel  jettisoning,  where  the  airstream  might 
not  possess  .suff  icient  energy  to  thoroughly  atomize  all  of  the  : . 

i  During  the  same  period  another  study  was  conducted  in  England  by 
j ------  ---------Cross  and  Ficknett  (Reference  6).  In  this  study  a  small  jet 

aircraft  (a  British  "Buccaneer")  jettisoned  fuel  at  7.5  kilograms 
.  per  second  from  wing  tanks  while  flying  at  120  meters  per  second. 
The  aircraft  was  flown  less  than  15  meters  above  the  ground,  and 
the  fuel  spray  was  sampled  by  filter  papers  laid  on  the  ground. 

A  fluorescent  dye  was  added  to  the  fuel  so  that  the  droplet  spots 
on  the  filter  paper  could  be  visualized.  The  fuel  used  was  Avtur, 
a  British  commercial  fuel  similar  to  JP-d  and  much  less  volatile 
S  than  JP-4.  A  typical  droplet  size  distribution  observed  for  one 

of  their  flights  is  shown  in  Figure  6  ("Cross  and  Picknett  - 
|  experimental").  The  mass  median  diameter  was  240  microns  with 

-  droplets  ranging  from  40  to  400  microns.  Cross  and  Picknett  con- 

?  sldered  this  distribution  to  be  a  lower  bound  for  fuel  jettisoning, 

E: 

f-  since  the  fuel  was  jettisoned  perpendicular  to  the  airstream.  As 

an  upper  bound,  they  presented  a  composite  distribution  based  on 
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previous  studies  for  other  aircraft  with  rearward  projecting 
jettison-tubes.  This  distribution  is  also  shown  in  Pigure  6 
("Cross  and  Picknett  -  Composite").  The  mass  median  diameter  in 
this  case  is  270  microns,  with  some  droplets  as  large  a  650 
microns.  ~ 

Although  the  aircraft-generated  droplet  size  distributions 
measured  by  Cross  and  Picknett  were  more  likely  to  be  represen¬ 
tative  of  aircraft  fuel  jettisoning  than  the  wind  tunnel  data  from 
AEDC,  there  was  still  some  concern  about  the  wide  disparity  of  the 
two  results.  Also,  the  Cross  and  Picknett  study  suffered  from  two 
shortcomings:  first,  the  technique  used  to  visualize  the  droplets 
could  not  detect  droplets  less  than  40  microns  in  diameter;  and 
second,  in  all  of  the  trials,  no  more  than  55  percent  of  the  jet¬ 
tisoned  fuel  could  ever  be  accounted  for.  In  an  attempt,  to 
resolve  the  discrepant  results,  Dawbarn  (Reference  4)  suggested 
that  the  actual  fuel  droplet  distribution  was  somewhere  between 
the  two  extremes.  His  proposed  intermediate  distribution  is  shown 
in  Pigure  6  as  "AEDC  -  Suggested  Compromise."  Due  to  the  uncer¬ 
tainty  still  surrounding  the  size  distribution  question,  the  deci¬ 
sion  was  made  to  perform  another  study  in  which  an  aircraft  would 
be  used  to  Jettison  the  fuel.  The  distinguishing  feature  of  this 
study  would  be  in-flight  sampling  by  a  second  aircraft. 

3.2  DESCRIPTION  OP  THE  EXPERIMENTAL  PROGRAM 

The • in-flight  sampling  program  was  conducted  jointly  by  the  Air 
Force  Engineering  and  Services  Center  and  the  Air  Force  Geophysics 


Laboratory,  at  Edwards  Air  Force  Base,  California,  in  December 
1976.  A  sampling  aircraft  equipped  with  optical  droplet-sizing 
instrumentation  and  a  hydrocarbon  vapor  analyzer  wa3  flown 
directly  into  the  cloud  of  fuel  jettisoned  by  an  Air  Force  KC-135 
tanker  aircraft.  Direct  measurements  of  the  fuel  droplet  size 
distribution  and  density  were  obtained  and  average  droplet  fall 
rates  were  estimated.  Flights  were  conducted  at  three  altitudes 
(approximately  1500,  3600,  and  6000  meters  above  mean  sea  level) 
to  determine  whether  the  fuel  droplet  size  distribution  was  a 
function  of  release  height.  Ground-level  sampling  for  hydrocar¬ 
bons  was  also  performed  during  the  low-level  flights  (which  were 
roughly  750  meters  above  the  ground)  to  determine  whether  any  of 
the  jettisoned  fuel  reached  the  ground.  In  all  cases,  typical 
KC-135  jettisoning  practices  were  followed.  Table  2  lists  the 
conditions  under  which  the  fuel  was  jettisoned  in  this  program. 


TABLE  2.  KC-135  FUEL  JETTISONING  PARAMETERS 


Fuel : 

Dump-nozzle  diameter: 
Aircraft  airspeed: 
Jettison  rate: 

Duration  of  dump: 


JP--4 

10.1  centimeters 

170  meters  per  second 

56  kilograms  per  second 
or  330  grams  per  meter 

15-30  seconds 


The  KC-135  jettisons  fuel  through  the  boom  normally  used  to  refuel 
other  aircraft.  Usually  the  boom  is  lowered  for  jettisoning  as 
well  as  refueling,  but  in  some  emergency  situations  fuel  might  be 


Jettisoned  with  the  boom  stowed.  Therefore,  the  experimental 
program  included  jettisoning  in  both  configurations. 

3.3  IN-FLIGHT  SAMPLING  RESULTS 

The  in-flight  sampling  was  carried  out  by  Meteorology  Research, 
Incorporated  (MRI),  using  a  specially  instrumented  Piper  Navajo 
aircraft.  The  droplet  sizing  instrumentation  consisted  of  three 
particle  spectrometers  manufactured  by  Particle  Measuring  Systems, 
Incorporated.  These  spectrometers,  which  size  particles  as  they 
intercept  a  beam  of  laser  light,  cover  the  size  ranges  2-30, 
20-300,  and  300-4500  microns.  A  foil  impactor  was  also  used  on 
some  flights  to  provide  a  larger  cross-section  for  detecting  drop 
lets  greater  than  250  microns  in  diameter.  Other  Instrumentation 
Included  a  hydrocarbon  analyzer  and  an  integrating  nephelometer  to 
detect  vapor  and  fine  particles,  respectively.  A  two-volume 
report  describes  the  MRI  sampling  effort  and  presents  in  detail 
the  data  obtained  (Reference  7)* 

After  the  tanker  performed  each  fuel  dump,  the  sampling  aircraft 
was  vectored  in  for  perpendicular  passes  through  the  fuel  dump 
plume.  In  the  later  sorties,  better  sampling  results  were 
obtained  by  flying  the  sampling  aircraft  directly  down  the  length 
of  the  fuel  dump.  Plying  at  80  meters  per  second,  the  aircraft 
was  then  able  to  obtain  continuous  particle  counts  for  as  long  as 
30  to  45  seconds  on  a  single  pass.  The  fuel  dump  could  be  seen 
initially  as  a  long,  white  cloud  resembling  a  contrail.  For 
favorable  sun  angles  the  fuel  cloud  was  still  visible  more  than  10 


minutes  after  the  dump,  but  usually  the  sampling  aircraft  could  no 
longer  follow  the  cloud  visually  after  the  first  two  to  four 
minutes.  Real-time  radar  tracking  of  the  tanker’s  position  during 
the  fuel  dump  and  of  the  location  of  the  sampling  aircraft  during 
successful  passes  through  the  plume  was  used  to  direct  the  search 
after  visual  sighting  was  no  longer  possible.  Unexpectedly, 
n-ur  drizzle-sized  droplets  (around  100  microns)  could  be  detected  even 
after  the  fuel  cloud  was  ro  longer  visible.  In  fact  "splash  on 
the  windshield"  of  the  Navajo  became  the  most  definitive  mark  of 
entry  into  the  fuel  cloud. 

Figure  7  shows  a  typical  plume  produced  by  fuel  jettisoning. 

During  the  first  minute  the  plume  quickly  expands  to  nearly  100 
meters  in  diameter  and  takes  on  a  markedly  colled  appearance  under 
the  influences  of  the  aircraft  wake  vortex.  This  long  spiral  tube 
maintains  Its  integrity  for  several  minutes,  but  eventually  begins 
to  fold  on  itself  as  the  plume  dissipates.  The  position  of  the 
refueling  boom  greatly  affects  the  initial  growth  of  the  fuel  dump 
plume.  When  the  boom  is  up  in  the  stowed  position,  the  fuel  Is 
immediately  engulfed  in  one  of  the  expanding  wing  tip  vortices, 
leading  to  rapid  expansion  and  coiling.  However,  with  the  boom 
lowered  the  plume  Initially  lies  below  the  aircraft  wake  and, 
therefore,  expands  more  slowly  with  little  sign  of  coiling. 
Finally,  after  approximately  60  seconds  the  plume  does  become 
entrained  in  one  of  the  wake  vortices  and  from  then  on  is  iden¬ 
tical  to  the  boom-up  plume.  The  pilot  of  the  Navajo  reported  that 
visually  there  appeared  to  be  greater  vertical  dispersion  of  the 
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plume  and  possible  "rainout"  of  larger  droplets  with  the  boom 
down,  but  no  sampling  or  photographic  evidence  of  this  effect 
could  be  obtained.  Sampling  passes  below  the  visible  fuel  clouds 
were  performed  to  search  for  any  rainout  of  larger  droplets,  but 
no  such  phenomenon  was  ever  detected. 

Overall  planning,  direction  and  analysis  of  tne  MRI  sampling 
effort  was  performed  by  the  Air  Force  Geophysics  Laboratory 
(AFGL).  Using  a  statistical  analysis,  AFGL  was  able  to  distin¬ 
guish  actual  droplet  encounters  from  fluctuations  in  the  atmos¬ 
pheric  particulate  background.  The  droplet  size  distribution 
—  observed  in  these  encounters  was  then  corrected  for  evaporation 
(using  the  droplet  evaporation  model  described  in  Section  II  of 
this  report)  to  obtain  an  estimate  of  the  "original"  droplet  size 
distribution.  This  correction  was  necessary  to  allow  for  the 
length  of  time  the  fuel  droplets  were  free  to  evaporate  before  the 
sampling  aircraft  entered  the  fuel  dump  plume.  Due  to  flying 
safety  considerations ,  the  first  pass  by  the  sampling  aircraft 
through  the  fuel  dump  usually  occurred  as  much  as  90  seconds  after 
the  tanker  had  passed.  Based  on  the  evaporation  model,  a  typical 
fuel  droplet  had  lost  more  than  80  percent  of  its  original  mass  by 
this  time  and  had  been  reduced  to  about  half  of  its  original  size. 

Composite  experimental  and  corrected  droplet  size  distributions 
based  on  all  of  the  sampling  flights  are  shown  in  Figure  8.  (The 
curve  entitled  "This  Report"  will  be  discussed  later.)  No  clear 
distinction  was  found  between  results  for  different  altitudes  and 
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boom  configurations;  therefore,  all  of  the  encounters  were  averaged 
together  to  obtain  these  composites.  Individual  sampling  results 
for  each  fuel  dump  plume  encounter  are  included  in  the  AFGL  report 
(Reference  8).  A  full  range  of  droplets  from  2-600  microns  were 
observed,  with  a  composite  mass  median  diameter  of  165  microns. 
After  correcting  for  evaporation,  the  mass  median  diameter  becomes 
265  microns,  with  a  significant  number  of  droplets  as  large  as  500 
microns. 

Ground-level  photography  was  used  by  AFGL  to  estimate  the  dimen¬ 
sions  of  the  fuel  dump  plume  during  the  sampling  operation.  Their 
analysis  indicates  that  the  fuel  dump  plumes  were  typically  100 
meters  in  diameter  after  90  seconds.  For  a  jettisoning  rate  of 
330  grama  per  meter,  the  nominal  density  of  fuel  dumped  would  then 
be  0.042  grams  per  cubic  meter.  Correcting  the  observed  fuel 
droplet  densities  for  evaporation,  AFGL  was  able  to  account  for 
only  0.010  grams  per  cubic  meter,  or  about  25  percent  of  the  known 
mass.  Similarly,  assuming  the  fuel  droplets  were  80  percent  evap¬ 
orated  at  the  time  the  plume  was  sampled,  the  measured  JP-4  vapor 
density  should  have  been  80  percent  of  the  nominal  fuel  density, 
or  0.034  grams  per  cubic  meter  (about  6.4  part  per  million).  The 
highest  concentration  recorded  by  the  hydrocarbon  analyzer  on  the 
Navajo  was  only  4  parts  per  million,  and  average  concentrations 
for  passes  down  the  length  of  the  plume  were  closer  to  2  parts  per 
million.  Since  both  of  these  discrepancies  diminished  confidence 
in  the  validity  of  the  sampling  results,  they  will  be  examined 
further  in  the  following  sections. 


3.3.1  ANALYSIS  OF  THE  DROPLET  MEASUREMENTS 


During  the  first  experimental  sortie,  only  perpendicular  passes 
through  the  fuel  dump  plume  were  performed.  However,  it  soon 
became  apparent  that  these  perpendicular  encounters  were  too  short 
to  provide  significant  particle  counts.  Thereafter,  parallel 
passes  down  the  length  of  the  fuel  dump  plume  were  also  carried 
out.  The  parallel  passes  are  preferred  not  only  because  they 
encountered  a  larger  number  of  particles,  but  also  because  they 
would  tend  to  average  out  any  local  fluctuations  in  the  fuel  dump 
plume  density.  The  longest  continuous  encounter  with  each  fuel 
dump  plume  is  listed  in  Table  3.  Fuel  dumps  not  listed  did  not 
include  any  long  sampling  passes. 

Two  of  the  fuel  dump  encounters  listed  in  Table  3  were  not  used 
In  the  mass  balance  calculation  that  follows.  In  fuel  dump  3/5, 
the  plume  rose  above  the  initial  jettisoning  height  under  the 
influence  of  a  high  level  Inversion  with  a  vertical  wind  com¬ 
ponent.  Apparently,  in  pushing  the  plume  upward,  the  inversion 
also  compressed  it  vertically,  because  the  observed  fuel  droplet 
density  was  considerably  greater  than  that  observed  in  any  of  the 
other  fuel  dumps.  In  fuel  dump  3/6,  the  only  lengthy  encounter 
with  the  plume  recorded  an  anomalous  droplet  distribtuion.  All  of 
the  observed  mass  was  in  very  large  droplets  (between  230  and  300 
microns).  Since  these  results  are  contradictory  to  the  obser¬ 
vations  of  all  the  other  fuel  dumps,  they  were  rejected. 

The  mass  densities  listed  in  Table  3  are  those  observed  by  the 
Cloud  Particle  Spectrometer  (CPS)  only.  This  instrument  sizes 
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droplets  with  diameters  ranging  from  10  to  310  microns  by  sorting 
the  droplets  into  15  channels,  each  of  which  includes  a  particle 
diameter  range  of  20  microns.  As  can  be  seen  from  the  experimen¬ 
tal  curve  in  Figure  8,  practically  all  of  the  observed  droplet 
mass  was  within  the  range  of  the  CPS.  Eliminating  consideration 
of  the  smaj.1  (2-20  micron)  and  large  (300-4500  micron)  particle 
instruments  simplifies  the  mass  balance  calculations  without 
significantly  changing  the  result. 

—-Calculation  of  the  average  intial  droplet  mass  density  for  the 
six  fuel  dumps  at  1500  meters  is  shown  in  Table  4.  The  average 
mass  density  observed  in  each  CPS  channel,  taken  directly  from 
Reference  7>  was  multiplied  by  a  correction  factor.  This  factor 
was  necessary  because  in  converting  from  cot^eotel  (raw  minus 
background)  droplet  counts  to  mass  of  JP-4,  MRI  used  the  wrong 
density  for  JP-4  (0.92  Instead  of  0.76).  Also,  their  conversion 
to  mass  neglected  the  effect  of  the  median  diameter  offset  which 
they  had  noted  for  the  first  five  channels.  These  discrepancies 
were  also  corrected  in  the  AFQL  analysis  (cf.  Table  2  in 
Reference  8  versus  Table  2  in  Reference  7).  The  fuel  droplet 
evaporation  model  was  used  to  determine  the  initial  droplet 
diameter  which  would  evaporate  to  each  channel's  median  diameter 
in  the  time  available  and  at  the  observed  ambient  temperature. 
Using  the  percent  mass  remaining  calculated  by  the  model,  the 
observed  mass  density  for  each  channel  was  then  scaled  up  to 
obtain  the  initial  droplet  mass  densities. 
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The  resulting  mass  balance  at  all  three  altitudes  is  shown  in 
Table  5*  The  ratio  of  the  calculated  initial  mass  density  to  the 
nominal  density,  based  on  a  plume  diameter  of  100  meters,  ranges 
from  0.6  to  1.8,  with  an  average  of  1.1.  The  composite  initial 
^distribution  for  all  ten  encounters  is  shown  in  Figure  8  ("This 
Report  -  Corrected").  This  distribution  differs  from  the  Initial 
distribution  derived  by  AFGL  only  in  the  smaller  droplet  diameter 
range.  Where  they  differ,  the  AFGL  distribution  is  probably  the 
more  accurate  for  two  reasons:  first,  the  AFGL  distribution  is 
based  on  data  obtained  by  all  three  spectrometers,  not  just  the 
CPS;  and  second,  AFGL  used  a  more  sophisticated  method  of  back¬ 
ground  subtraction  to  obtain  the  corrected  droplet  counts  than 
that  used  by  MR I  in. Reference  7. 

The  total  initial  density  derived  by  AFGL  was  only  0.010  grams  per 
cubic  meter,  compared  to  an  average  of  0.046  derived  here.  How¬ 
ever,  the  curves  of  initial  droplet  density  versus  diameter  in  the 
AFGL  report  (Reference  8)  are  consistent  with  the  results  obtained 
here.  For  example,  the  initial  droplet  density  curve  in  the  AFGL 
report  for  the  fuel  dumps  at  1500  meters  rises  from  1.5  x  10 
grams  per  cubic  meter  at  100  microns  to  a  maximum  of  8.3  x  10 
grams  per  cubic  meter  at  350  microns,  and  then  falls  to  around 
1  x  10  grams  per  cubic  meter  at  500  microns.  From  Table  4,  the 

_-p 

corresponding  densities  are  1.7  x  10  grams  per  cubic  meter  at 
100  microns,  4.7  x  10“^  grams  per  cubic  meter  at  340  microns,  and 
1.41  x  10  grams  per  cubic  meter  at  500  microns.  The  discrep¬ 
ancy,  then,  only  appears  when  these  droplet  density  curves  are 
Integrated  to  obtain  the  total  liquid  density. 


TABLE  5.  SUMMARY  OF  DROPLET  MASS  DENSITY  CALCULATIONS 


As  can  be  seen  from  Table  4,  in  the  process  of  correcting  for 
evaporation  the  significance  of  the  actual  instrument  channel 
width  is  lost,  and  the  question  of  what  "bin  size"  to  use  in  the 
integration  arises.  The  analysis  performed  here  was  able  to 
simply  add  the  results  for  each  channel  of  the  CPS.  Since  the 
AFGL  analysis  was  based  on  three  overlapping  spectrometers  with 
different  channel  widths,  it  was  necessary  to  perform  the  integra¬ 
tion  analytically.  To  accomplish  this,  expressions  for  the 
droplet  density  as  a  function  of  diameter  (D)  were  fit  to  the 
curves.  Based  on  discussions  with  the  AFGL  project  officer,  these 
expressions  were  then  divided  by  D  and  integrated,  yielding  the 
total  of  0.010  grams  per  cubic  meter.  The  usual  procedure  would 
be  to  divide  the  density  function  by  the  bln  size  (if  it  were 
known)  rather  than  by  D.  If  the  density  functions  from  Reference 
b  are  integrated  and  divided  by  an  "effective"  bin  size  of  30 
microns  (from  Inspection  of  Table  4),  the  total  initial  liquid 
density  becomes  O.Obb  grams  per  cubic  meter.  This  alternative 
calculation  is  somewhat  arbitrary,  but  the  resultant  density  is 
more  in  line  with  the  results  in  Table  5,  and  the  proper  way  to 
arrive  at  an  effective  bin  size  for  this  case  is  not  known.  At 
any  rate,  the  observed  droplet  mass  densities,  based  on  CPS  data, 
appear  to  be  within  expectations. 

As  shown  in  Figure  9,  the  initial  fuel  droplet  size  distribution 
obtained  for  the  KC-135  in  this  study  agrees  very  well  with  the 
distributions  reported  by  Cross  and  Picknett  (Reference  6).  The 
smaller  droplet  size  distribution  found  by  AEDC  (Reference  5)  is 
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evidently  not  representative  of  full  scale  Jettisoning  by  aircraft. 
The  agreement  with  Cross  and  Picknett 's  results  is  even  more 
-remarkable  if  the  possible  effect  of  evaporation  on  their  results 
is  considered.  As  mentioned  earlier.  Cross  and  Picknett  could 
never  account  for  more  than  55  percent  of  the  Jettisoned  fuel.  To 
investigate  the  extent  of  evaporation,  they  measured  the  concen¬ 
tration  of  dye  in  the  spots  on  the  filter  paper  and  compared  this 
with  the  amount  of  dye  added  to  the  fuel.  On  one  trial  it  appeared 
that  drops  initially  117  microns  in  diameter  had  evaporated  to  100 
microns  before  being  sampled,  but  they  were  unable  to  confirm  this 
-result  in  two  subsequent  trials.  They  therefore  concluded  that  no 
significant  evaporation  took  place.  However,  a  change  in  size 
from  117  to  100  microns  amounts  to  the  loss  of  40  percent  of  the 
droplet's  mass.  Therefore,  to  further  examine  the  possibility  of 
evaporation,  the  fuel  droplet  evaporation  model  was  run  for  the 
conditions  of  their  experiment. 

Since  Cross  and  Picknett' s  study  was  performed  in  August  in 
England,  an  ambient  temperature  of  25°C  was  assumed.  To  simulate 
the  relatively  non-volatile  Avtur  fuel  used  in  their  study,  a  com¬ 
position  based  on  analyses  of  JP-8  and  Jet  A-l  was  input  to  the 
model  instead  of  the  synthetic  JP-4  composition  shown  in  Table  1. 
The  model  was  then  run  starting  with  the  fuel  droplet  size  distri¬ 
bution  found  for  the  KC-135.  For  a  release  height  of  15  meters, 
the  model  predicts  that  approximately  40  percent  of  the  fuel  will 
have  evaporated  before  it  reaches  the  ground.  The  predicted 
droplet  distribution  at  the  ground  is  shown  in  Figure  9  ("KC-135 


Distribution  -  After  40  Percent  Evaporated”).  Although  the  mass 
median  diameter  of  the  distribution  is  not  significantly  altered 
by  the  evaporation,  the  relative  contribution  from  smaller 
droplets  is  greatly  reduced.  Thus  allowing  for  the  possibility  of 
evaporation  in  Cross  and  Picknett's  experiment  not  only  accounts 
for  their  missing  mass,  but  also  reduces  the  difference  between 
their  distribution  and  the  distribution  observed  in  this  study. 
Based  on  the  agreement  between  these  two  studies  it  is  postulated 
that  the  distribution  obtained  for  the  KC-135  aircraft  is  repre¬ 
sentative  of  most  other  aircraft  as  well. 

3.3.2  ANALYSIS  OP  VAPOR  MEASUREMENTS 

The  average  hydrocarbon  vapor  concentrations  measured  during  long 
passes  through  the  fuel  dumps  were  a  factor  of  two  to  three  below 
the  levels  calculated  assuming  80  percent  evaporation  of  the  fuel 
and  a  plume  diameter  of  100  meters.  Since  the  droplet  measure¬ 
ments  have  been  shown  to  be  consistent  with  both  of  these  assump¬ 
tions,  the  discrepancy  is  probably  due  to  some  deficiency  in  the 
vapor  measurements.  Therefore,  the  adequacy  of  the  technique  used 
was  re-examined. 

The  instrument  used  to  obtain  the  hydrocarbon  measurements  was  a 
Century  Portable  Organic  Vapor  Analyzer  (OVA)  Model  108.  This 
instrument  features  a  logarithmic  scale  covering  the  range  from  1 
to  10,000  parts  per  million  and  a  response  time  under  two  seconds. 
For  use  in  this  study  an  in-line  carbon  filter  was  removed  from 
the  inlet  hose  and  the  inlet  hose  was  then  positioned  in  an 


ambient  air  inlet  on  the  skin  of  the  aircraft.  A  reverse  flow 
inlet  was  used  to  avoid  sampling  particulates  larger  than  60 
microns.  Also,  to  assure  continuous  operation  of  the  OVA’s  flame 
ionization  detector  at  the  altitudes  included  in  this  study,  it 
was  necessary  to  increase  the  hydrogen  supply  pressure  from  8 
pounds  per  square  inch  to  15  pounds  per  square  inch. 

Surprisingly,  this  modification  resulted  in  about  a  ten-fold 
Increase  in  response  to  methane  standards, 

gBefore  the  in-flight  sampling  program  was  performed,  the  OVA  was 
calibrated  in  an  altitude  chamber  at  each  of  the  planned  altitudes 
Aluminum  clyinders  containing  three  concentrations  of  methane  in 
dry  air  were  used:  10,  102.4,  and  1054  parts  per  million.  These 
standards  were  transferred  to  polyethlyene  bags  for  sampling. 

These  studies  indicated  a  reduction  in  response  of  approximately 
20  and  40  percent  in  going  from  sea-level  to  3600  and  6000  meters, 
respectively.  However,  actual  in-flight  calibration  with  the  same 
standards  indicated  much  more  severe  reductions  in  response:  73 
percent  at  3600  meters  and  96  percent  at  6000  meters.  The  con¬ 
tractor  assumed  this  difference  was  due  to  the  low  ambient  tem¬ 
peratures  encountered  during  the  in-flight  sampling.  Therefore 
this  "temperature  effect"  was  incorporated  into  the  OVA  calibra¬ 
tion  equation. 

The  relative  response  of  the  OVA  to  JP-4  vapor  was  determined  by 
injecting  a  known  amount  of  JP-4  liquid  into  a  plastic  bag  and 
allowing  the  J-P4  to  evaporate.  In  calculating  the  concentration 


thus  obtained,  the  contractor  unfortunately  used  an  incorrect 
value  for  the  density  and  mean  molecular  weight  of  JP-*!.  Instead 
of  0.92  grams  per  milliliter  and  170  grams  per  mole,  the  contrac¬ 
tor  should  have  used  0.76  grams  per  milliliter  and  127  grams  per 
mole.  However,  the  effects  of  these  two  errors  partially  cancel 
and  the  overall  effect  is  to  raise  the  concentrations  reported  in 
—  Reference  7  by  only  9  percent.  The  relative  response  of  the  OVA 
to  JP-4  vapor,  as  determined  by  the  above  method,  is  92  percent  of 
the  response  to  methane. 

After  the  sampling  program,  when  it  was  found  that  observed 
hydrocarbon  vapor  levels  had  been  less  than  expected,  possible 
explanations  for  the  discrepancy  were  studied.  Although  the 
actual  instrument  used  in  the  in-flight  sampling  was  no  longer 
"available,  an  identical  OVA  Model  108  was  obtained  and  tested. 

A  relatively  flat  response  to  hydrocarbons  of  varying  molecular 
weight  was  verified,  eliminating  any  concern  about  differences 
between  the  composition  of  the  vapor  sampled  in  flight  and  that 
produced  during  the  calibration.  The  increase  in  response 
reported  by  the  contractor  when  the  hydrogen  supply  pressure  was 
raised  was  also  reproduced.  However,  the  decrease  in  response 
with  temperature  suggested  by  the  contractor  to  explain  the  dif- 
f .  ce  between  the  In-flight  and  altitude  chamber  calibration 
results  was  not  found.  The  response  of  this  second  OVA  to  methane 
and  propane  standards  was  unaffected  by  lowering  the  ambient  tem¬ 
per'  re  from  23°C  to  2°C.  Of  course,  this  does  not  completely 
exclude  the  possibility  that  the  instrument  used  by  MRI  did  indeed 
possess  a  temperature  dependence. 
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The  in-flight  calibrations  were  performed  by  filling  plastic  bags 
from  the  standards  in  aluminum  cylinders.  The  OVA  inlet  hose  was 
then  removed  from  the  air  inlet  on  the  skin  of  the  aircraft  and 
introduced  into  the  bag.  If,  as  suggested  by  the  contractor,  the 
response  of  their  OVA  was  a  function  of  temperature,  then  the  in¬ 
flight  calibration  would  not  be  completely  accurate.  Although  the 
sampling  aircraft’ 3  heater  was  not  on,  the  warming  effect  of  the 
personnel  and  equipment  on  board  would  tend  to  keep  the  interior 
of  the  aircraft  well  above  ambient  temperatures,  particularly  at 
the  higher  altitudes  where  outside  temperatures  were  below  0°C. 

Thus  the  response  of  the  instrument  to  the  standards  inside  the 
aircraft  would  be  greater  than  its  response  to  the  JP-4  vapor  out¬ 
side  where  it  was  colder.  This  effect  would  then  explain  the 
failure  to  observe  the  anticipated  vapor  levels.  On  the  other 
hand,  if  the  response  of  the  OVA  was  not  affected  by  temperature, 
as  was  the  case  for  the  second  OVA  studied  later,  then  there  is  no 
explanation  why  the  response  of  the  OVA  during  in-flight  calibra¬ 
tion  was  so  much  lower  than  during  the  altitude  chamber  runs.  In 
this  case  it  could  only  be  assumed  that  some  unidentified  instru¬ 
ment  malfunction  occurred,  compromising  the  in-flight  measurements. 

3.3.3  EFFECT  OF  IGNORING  THE  WAKE  IN  THE  MODEL 

During  the  first  few  minutes  after  being  jettisoned,  fuel  droplets 
are  subject  to  the  complex  forces  present  in  the  aircraft  wake. 
Therefore,  the  droplets  are  not  in  "free-fall"  until  the  wake 
dissipates.  The  fuel  droplet  evaporation  and  free-fall  model. 
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however,  assumes  that  from  the  start  gravitational  settling  is  the 
only  force  acting  on  the  droplets.  Nevertheless,  this  model  was 
used  to  estimate  the  evaporation  which  occurred  during  the  first 
few  minutes  before  the  droplets  were  sampled,  in  order  to  arrive 
at  an  initial  droplet  size  distribution.  It  is  reasonable  then  to 
question  the  propriety  of  ignoring  the  effect  of  the  wake  in  the 
model.  This  question  will  be  considered  from  two  viewpoints: 
first,  the  empirical  evidence  for  the  sufficiency  of  the  model  is 
given;  and  second,  the  theoretical  basis  for  this  sufficiency  is 
suggested. 


Prom  the  standpoint  of  using  the  model  to  predict  the  later 
history  of  the  fuel  droplets,  the  accuracy  of  the  calculated  ini¬ 
tial  distribution  is  of  secondary  importance.  Since  the  model 
that  was  used  to  work  backwards  from  the  observed  distribution  at 
90  seconds  to  an  Initial  distribution  is  the  same  model  used  to 
obtain  the  later  predictions,  the  purpose  of  obtaining  an  initial 
distribution  is  only  to  assure  that  the  model  passes  through  the 
observed  distribution  and  mass  remaining  at  the  obervatlon  time. 
After  this  point  the  droplets  are  essentially  in  free  fall  and  the 
model  assumptions  are  valid.  The  calculated  fuel  droplet  distri¬ 
butions  at  90  seconds  (the  typical  observation  time)  for  the  three 
altitudes  Included  in  this  study  are  compared  in  Figure  10  with 
the  experimentally  observed  distributions.  The  calculated  initial 
distribution  is  also  shown  for  contrast.  At  all  three  altitudes 
the  agreement  is  considered  satisfactory.  The  percentage  of  mass 
remaining  at  90  seconds  calculated  by  the  droplet  evaporation 


model  for  the  1500,  3600,  and  6000  meter  cases  is  12,  15,  and  30 
percent,  respectively.  The  corresponding  observed  values,  baaed 
on  a  nominal  dump  rate  of  330  grams  per  meter  and  a  plume  diameter 
of  100  meters  at  90  seconds,  were  14,  26,  and  17  percent.  Once 
again,  this  rough  agreement  is  considered  satisfactory,  given  the 
variability  of  the  experimental  results. 

In  order  to  explain  the  apparent  similarity  of  droplet  evaporation 
in  the  wake  and  in  free-fall,  it  is  necessary  to  consider  the  for- 
ces  involved  in  each  case.  As  discussed  in  Appendix  A,  the  'rate 
of  evaporation  of  a  droplet  is  the  product  of  the  vapor  pressure 
of  its  components,  the  surface  area  available  for  evaporation,  and 
the  mass  transfer  coefficient.  The  mass  transfer  coefficient 
depends  in  part  on  the  droplet’s  terminal  velocity,  which  is 
determined  by  the  external  forces  acting  on  the  droplet.  Whether 
the  droplet  is  experiencing  the  wake  vortex  or  just  gravitational 
settling,  an  equilibrium  is  reached  between  the  forces  accelerating 
the  droplet  and  the  drag  resisting  its  motion.  This  equilibrium 
determines  the  droplet's  terminal  velocity  and  thus,  indirectly, 
its  evaporation  rate.  During  free-fall  the  acceleration  is  simply 
that  due  to  gravity,  or  9,8  meters  per  square  second.  For  a 

p 

stable  orbit  in  the  wake,  the  centripetal  acceleration  is  V  /r, 
where  V  is  the  angular  velocity  in  the  wake,  and  r  is  the  orbit 
radius.  This  stable  orbit  is  maintained  by  the  strong  influx  of 
air  being  entrained  into  the  growing  vortex. 
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TABLE  6.  ESTIMATION  OP  TYPICAL  DROPLET  ACCELERATION  IN  WAKES 


Orbit 

Angular 

Centripetal 

Radius 

Velocity 

Acceleration 

(m/s2) 

Aircraft 

(meters) 

( m/  s ) 

DC-6B 

1 

8 

64 

(wake  age  ■  25 

sec)  4 

4 

4 

6 

3 

1.5 

8 

3 

1.1 

12 

2 

0.3 

0-141 

1 

15 

225 

(wake  age  *  10 

sec)  2 

11 

60 

4 

7 

12 

6 

6 

6 

8 

5 

3 

"w  : :  - 

12 

3 

0.8 

C-141 

2 

5 

12 

(older  wake) 

14 

3 

0.6 

Table  6  shows  typical  values  of  acceleration  that  a  droplet  can 
experience  In  an  aircraft  wake  (Reference  9).  The  aircraft  shown 
are  roughly  the  same  size  and  weight  as  a  KC-135>  and  therefore 
the  wake  energies  should  be  similar.  Although  the  range  of  accel¬ 
eration  appears  large,  the  net  effect  on  the  droplet  evaporation 
rate  is  small.  If  the  droplet  evaporation  calculations  in 
Appendix  A  are  carried  out  using  the  centripetal  accelerations 
from  Table  6  instead  of  the  gravitational  acceleration  (g),  the 
resulting  evaporation  rates  range  from  0.6  to  2.1  times  the  rate 
for  free-fall,  with  typical  values  lying  near  unity.  In  other 
words,  the  average  environment  experienced  by  fuel  droplets  while 
in  the  aircraft  wake  does  not  significantly  alter  their  evapora¬ 
tion  rate  from  that  which  occurs  during  simple  gravitational 
settling. 
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Ignoring  the  aircraft  wake  does  substantially  underestimate  the 
initial  fall  rate  of  the  droplets,  however.  From  windshield 
splash  encounters  and  radar  records  the  initial  fall  rate  of  the 
fuel  droplets  was  around  150  centimeters  per  second,  but  after 
approximately  100  seconds  the  observed  fall  rate  had  reduced  to 
around  35-50  centimeters  per  second  (Reference  8).  The  faster 
Initial  fall  rate  is  due  to  the  downward  force  of  the  aircraft 
wake.  Thereafter,  the  droplets  are  in  free-fall  and  the  observed 
fall  rate  is  consistent  with  the  predictions  of  the  model  for 
droplets  on  the  order  of  120  to  160  microns  in  diameter.  Thus, 
ignoring  the  wake  effect  only  underestimates  initial  droplet 
fall-rates,  and  the  overall  error  of  about  100  meters  is  not 
significant  for  most  uses  of  the  model. 

3.4  GROUND-LEVEL  SAMPLING 

In  parallel  with  the  In-flight  sampling  program,  a  separate  pro¬ 
gram  of  ground-level  sampling  was  also  performed.  Whereas  the 
purpose  of  the  In-flight  sampling  was  to  measure  the  Initial  fuel 
droplet  size  distribution,  the  ground-level  sampling  was  designed 
to  determine  whether  any  of  the  jettisoned  fuel  reached  the 
ground  in  significant  concentrations,  either  as  vapor  or  liquid. 
For  this  purpose  18  sampling  stations  were  set  up  in  a  two-mile 
square  grid,  as  shown  in  Figure  11.  Several  of  the  low-altitude 
fuel  dumping  sorties  for  the  in-flight  sampling  were  positioned  so 
that  the  fuel  would  reach  the  ground  in  the  area  of  the  sampling 
grid.  Air  samples  were  obtained  with  small  diaphram  pumps,  using 
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charcoal  adsorption  tubes  to  collect  any  hydrocarbons  present. 

The  hydrocarbons  were  later  desorbed  from  the  charcoal  in  the 
laboratory  and  analyzed  by  gas  chromatography. 

The  positioning  of  the  fuel  dumps  was  excellent,  and  in  all  cases 
personnel  on  the  ground  at  the  sampling  grid  were  able  to  smell 
JP-4  vapor  approximately  13  minutes  after  release.  However,  none 
of  the  ground  observers  were  ever  able  to  detect  any  sign  of 
liquid  fuel  at  ground-level.  These  observations,  which  are  for 
fuel  dumped  approximately  750  meters  above  the  ground  at  tempera¬ 
tures  around  11°C,  substantiate  the  prediction  of  the  fuel  droplet 
evaporation  and  free-fall  model  that  under  such  conditions  essen¬ 
tially  all  of  the  fuel  will  evaporate  before  it  reaches  the 
ground. 

3.4.1  PREDICTED  CONCENTRATION  AND  LOCATION  OP  THE  FUEL  AT 
THE  GROUND 

Four  ground-sampling  runs  were  performed  in  all.  During  the 
second  in-flight  sampling  sortie  on  6  December,  ground  samples 
were  obtained  for  fuel  dumps  2/1  and  2/3  (using  the  notation  from 
Reference  8).  On  10  December,  fuel  dumps  3/2  and  3/3  were  sampled 
together  and  a  final  dump,  which  will  be  referred  to  as  3/7,  was 
performed  (with  the  boom  down)  for  the  ground-level  sampling 
alone.  The  winds  and  temperatures  at  the  time  of  the  sampling  are 
shown  in  Table  7. 


TABLE  7.  ATMOSPHERIC  CONDITIONS  DURING  GROUND-SAMPLING 

STUDIES 


Date  and 

Local  Time 

Altitude 
(meters,  MSL) 

Temperature 

(C) 

Wind 

Direction 

(degrees) 

Wind 

Spec 

(m/s 

6  Dec,  1430 

1  :  . 

1500 

7-5 

62 

7.3 

1200 

8.4 

65 

5.0 

900 

12.7 

53 

4.0 

725* 

15.4 

45 

3.5 

10  Dec,  1300 

-1500 

4.5 

61 

8.5 

1200 

7.1 

72 

7.0 

900 

11.1 

72 

4.5 

725* 

14.3 

70 

3.0 

*local  ground-level 


Based  on  a  Gaussian  dispersion  model,  the  maximum  concentration 
of  fuel  at  ground-level  can  be  expected  to  occur  when  the  vertical 
standard  deviation  (oz)  of  the  fuel  cloud  is  equal  to  0.707  H, 
where  H  is  the  height  of  the  release  above  the  ground  (Reference 
10).  In  this  case  H  is  775  meters,  making  the  desired  °z  548 
meters.  For  a  sunny  winter  day  in  Southern  California  with  surface 
winds  around  3  meters  per  second,  the  atmospheric  stability  cate¬ 
gory  is  probably  B  or  C.  From  Reference  10,  for  stability  B-C, 
az  will  equal  548  meters  when  the  plume  has  been  blown  approxima¬ 
tely  7  kilometers  downwind,  and  at  this  point  the  width  of  the 
central  plume  at  the  ground  will  be  at  least  1500  meters.  That 
is,  the  maximum  concentration  of  fuel  at  ground  level  will  occur  7 
kilometers  downwind  of  the  aerial  dump  location. 

Figure  11  shows  the  aerial  location  of  each  of  the  fuel  dumps  and 
the  predicted  location  of  the  maximum  ground-level  concentrations 
for  each.  The  locations  of  the  ground  imprints  were  determined 
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by  integrating  the  effects  of  the  wind  direction  and  speed  from 
the  release  height  to  the  ground.  The  maximum  ground  Imprints 
for  the  fuel  dumps  performed  on  10  December  cross  the  northwest 
corner  of  the  sampling  grid.  The  maximum  ground  concentrations 
for  the  fuel  dumps  performed  on  6  December  fall  beyond  the  grid. 
However,  initial  contact  of  the  fuel  cloud  with  the  ground  would 


occur  well  before  the  point  where  the  maximum  concentration  was 
reached.  Therefore,  some  of  the  fuel  could  be  expected  to  reach 
the  sampling  grid,  as  is  supported  by  the  ability  of  personnel  in 
the  sampling  area  to  smell  JP-4  after  each  dump. 


Based  on  the  fuel  droplet  evaporation  and  free-fall  model,  the 
fuel  will  have  been  more  than  99  percent  evaporated  before  it 
reaches  the  ground.  Nevertheless,  the  droplets  that  still  remain 
are  predicted  to  reach  the  ground  after  about  20  minutes.  For  the 
wind  profiles  shown  in  Table  7,  in  20  minutes  the  droplets  will 
be  blown  roughly  7  kilometers  downwind.  Thus  the  maximum  ground 
imprints  shown  in  Figure  11  are  indicative  of  droplet  impact  as 
well.  While  the  fuel  dumps  on  10  December  are  better  positioned 
to  look  for  the  predicted  maximum  ground  Impact,  those  on 
6  December  serve  as  a  check  for  anomalously  high  ground  contamina¬ 
tion.  If  the  fuel  were  to  reach  the  ground  sooner  than  expected, 
either  from  slower  evaporation  leading  to  faster  fall  rates  or 
from  an  unexpected  mechanism  mixing  the  fuel  to  the  ground 
prematurely,  much  greater  ground  impact  would  result.  However,  in 
this  event,  the  early  ground-fall  could  be  detected  in  dumps  2/1 
and  2/3* 
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The  highest,  ground  concentrations  which  can  be  anticipated  for 
these  fuel  dumps  can  be  calculated  following  the  procedure  from 
Reference  10 : 

Xmax  =  Q/iroyazUe 
Where  Q  =  the  emission  rate 

=  56  kilograms  per  second 
oy  =  the  horizontal  standard  deviation 
=  730  meters  at  7  kilometers 
az  =  the  vertical  standard  deviation 
=  5^8  meters  at  7  kilometers 
U  =  the  windspeed  (here  the  aircraft  speed) 

=  175  meters  per  second 

The  maximum  concentration  calculated  in  this  way  is  0.11  milli¬ 
grams  per  cubic  meter.,  which  is  equivalent  to  0.021  parts  per 
million  of  JP-4,  or  0.19  parts-per-million-carbon  (ppmC).  The 
duration  of  this  concentration  at  a  sampling  point  depends  on  the 
length  of  the  fuel  plume  and  the  rate  at  which  the  wind  sweeps  the 
plume  past  the  sampling  point.  The  longest  Imprint  of  a  fuel 
plume  on  the  sampling  grid  is  about  3  kilometers.  For  a  surface 
wind  of  3  meters  per  second,  the  duration  would  be  roughly  17 
minutes.  The  shortest  sampling  times  used  in  this  study  were 
around  30  minutes.  Therefore,  the  highest  time-weighted  average 
concentration  that  could  be  obtained  with  the  charcoal  tubes  is 
predicted  to  be  on  the  order  of  0.1  ppmC. 

3.4.2  DESCRIPTION  OF  THE  ANALYTICAL  TECHNIQUE 
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Small  diaphragm  pumps  of  the  type  used  for  personnel  monitoring 
were  used  to  draw  air  samples  through  standard  charcoal  sampling 
tubes.  The  sampling  rate  of  each  pump  was  calibrated  to  1.0 
liters  per  minute  using  a  soap-bubble  burette.  The  pumps  were 
started  before,  or  within  ten  minutes  after,  the  fuel  was 
jettisoned,  and  were  allowed  to  run  until  approximately  30  minutes 
after  the  jettisoning.  On  10  December,  background  samples  were 
also  obtained  before  any  fuel  had  been  jettisoned.  - 

The  hydrocarbons  adsorbed  on  the  charcoal  were  analyzed  by  gas 
chromatography.  The  technique  used  to  transfer  the  hydrocarbons 
from  the  charcoal  to  the  gas  chromatograph  is  shown  in  Figure  12. 
The  charcoal  tube  was  flushed  with  inert  gas  and  was  heated  to 
approximately  200°C  in  a  heater  block  to  desorb  any  hydrocarbons 
present.  The  desorbed  hydrocarbons  were  trapped  in  a  sample  loop 
cooled  to  liquid  nitrogen  temperatures  before  the  inert  gas  stream 
was  vented  to  the  atmosphere.  To  analyze  the  sample,  the  gas  flow 
through  the  sample  loop  was  diverted  into  the  chromatograph,  while 
the  sample  loop  itself  was  heated  with  a  hot  air  gun  to  revolati¬ 
lize  the  trapped  hydrocarbons.  The  hydrocarbons  were  then 
separated  on  a  Porapak-Q  micro-column  programmed  from  50°  to  240°C 
at  5°C  per  minute.  The  flame  ionization  detector  response  was 
output  to  an  electronic  integrator  to  provide  quantification. 

Standards  containing  13  ppmC  of  JP-4  vapor  were  prepared  by 
injecting  0.5  microliters  of  liquid  JP-4  into  a  50  liter  bag  and 
allowing  the  fuel  to  vaporize.  A  typical  gas  chromatogram  is 
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Figure  12.  Diagram  of  Charcoal  Tube  Desorption  Technique 
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shown  in  Figure  13a.  The  three  major  peaks  have  retention  times 
of  approximately  60,  240,  and  860  seconds.  The  response  and 
standard  deviation  for  the  calibration  runs  are  shown  in  Table  8. 
Due  to  the  excessive  variability  of  the  peak  at  860  seconds,  this 
peak  was  not  used  to  quantify  the  experimental  results. 

TABLE  8.  CALIBRATION  OF  HYDROCARBON  ANALYSIS 


Total  Area: 

Major  Peaks 

60  seconds: 
240  seconds: 
860  seconds: 


Response 

( area/ppmC/li+.er ) 
110  +  12 


24  +  4 
31  +  9 
11  +  8 


The  results  for  the  four  ambient  backgrounds  collected  on  10 
December  appear  in  Table  9«  Three  of  these  samples  resembled 
Figure  13b.  Note  that  all  of  the  major  peaks  for  JP-4  are  also 
present  in  ambient  air.  The  background  taken  at  T4  had  several 
additional  peaks,  notably  ones  at  180,  210,  430,  and  730  seconds. 
Many  of  the  samples  collected  during  the  ground-sampling  on  both 
6  and  10  December  contained  these  additional  peaks,  and  they  were 
also  present  in  ambient  air  samples  collected  in  the  Los  Angeles 
basin  as  well  as  near  the  flightline  at  Tyndall  Air  Force  Base. 

A  typical  chromatogram  is  shown  in  Figure  13c.  It  appears  that 
these  additional  peaks  may  be  indicative  of  hydrocarbon  pollutants 
from  automobile  and  aircraft  exhaust.  Their  appearance  in  some  of 
the  ground  samples  on  6  and  10  December  could  be  due  to  the 
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Figure  13.  Typical  Gas-Chromatograms  of  Hydrocarbon  Samples 


exhaust  from  the  trucks  used  to  get  between  sampling  stations. 

The  background  in  the  Los  Angeles  basin  calculated  for  the  peaks 
of  interest  is  also  shown  in  Table  9*  Fortunately,  although  the 
presence  of  these  additional  hydrocarbons  does  of  course  increase 
the  total  hydrocarbon  concentration,  it  does  not  appear  to  affect 
the  peaks  at  60  and  240  seconds.  These  two  peaks  were  therefore 
used  to  analyze  the  ground-sampling  results. 


TABLE  9.  AMBIENT  BACKGROUND  SAMPLES 


Sampling 

Background 

Concentration 

(ppmC)  Based 

On : 

Location 

Total  Area 

Peak  60 

Peak  240 

Peak  860 

T1 

0.70 

0.70 

0.43 

1.07 

T2 

0.63 

0.42 

0.29 

2.15 

T4# 

3.46 

1.06 

0.37 

4.62 

T7 

1.19 

1.15 

0.50 

3.86 

Average 

1.50 

0.83 

0.40 

2.93 

Pasadena* 

t 

0.29 

0.53 

13-35 

^includes  additional  peaks  associated  with  automobile  exhaust, 
tintegrator  overflow 


3.4.3  ANALYSIS  OF  SAMPLING  RESULTS 


The  results  of  the  four  sampling  runs  are  summarized  in  Table  10. 
Only  on  the  last  run  were  concentrations  greater  than  the  back¬ 
ground  obtained.  From  the  negative  results  on  6  December  we  can 
conclude  that  there  is  no  evidence  for  the  premature  arrival  of 
the  fuel  at  the  ground  and  the  increased  ground  concentrations 
which  would  result.  The  fact  that  higher  concentrations  were 
observed  on  the  second  10  December  run  but  not  the  first  may  be 
due  to  the  difference  in  sampling  times.  While  all  the  other 
samples  and  backgrounds  were  obtained  over  30-50  minutes,  the 


samples  for  the  first  run  on  10  December  were  collected  for  90 
minutes  in  order  to  encompass  two  fuel  dumps.  However,  later 
studies  indicated  that  increasing  sampling  time  from  50  to  90 
minutes  decreased  the  apparent  concentrations  by  roughly  a  factor 
of  three.  This  decrease  is  probably  due  to  saturation  of  the 
charcoal  or  displacement  of  the  hydrocarbons  by  other  compounds  in 
ambient  air. 

TABLE  10.  SUMMARY  OP  GROUND-SAMPLING  RESULTS 


Background  Concentration  (ppmC)  Based  On: 

Total  Area  Peak  60  Peak  240 


6  December  - 
average : 

fuel  dump  2/1 

0.57 

0.18 

0.14 

maximum : 

1.02 

0.25 

0.24 

6  December  - 
average : 

fuel  dump  2/3 

0.35 

0.26 

0.37 

maximum : 

0.75 

0.26 

0.37 

10  December 
and  3/2 
average : 

-  fuel  dump  3/1 

0.62 

0.45 

0.19 

max imum : 

0.96 

0.77 

0.37 

10  December 
average : 

-  fuel  dump  3/7 

2.04 

1.77 

0.53 

maximum : 

4.79 

3.15 

1.09 

The  chromatogram  for  station  T4  obtained  in  the  final  run  is 
shown  in  Figure  13d.  The  quantitative  results  from  this  run  for 
all  of  the  stations  are  listed  Table  11.  Here  the  average 
background  levels  from  Table  9  have  been  subtracted  out.  Concen¬ 
trations  greater  than  the  variation  of  the  background  have  been 
underlined  for  emphasis.  Depending  on  the  peak  used,  concentra¬ 
tions  on  the  order  of  several  tenths  of  a  ppmC  or  even  as  high  as 
2  ppmC  were  seen  at  several  sampling  stations.  The  highest 
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Tabl®T1iuF  Jp“^  CONCENTRATIONS 

A1  THE  ^OUND  PROM  FUEL  DUMP  3/7  S 


Total  Hydrocarbon 

OOnn^nf  ho  _ 


station 

Background : 

Sample  minus 

~ ~ - (.ppmc) 

(ppmC 

Peak 

O  Based  on: 
-§2  Peak 

1*5  +  2. 0f 

background : 

0.8  + 

0.4f  0.4  + 

TO 

- - -  .  . . 

T2 

0.4 

- 

T3# 

1.0 

0.8 

T4* 

1.3 

1. 6 

0.1 

T6# 

2.  0 

0.4 

T7* 

1.3 

1,8 

0.3 

T8* 

|u> 

2.4 

0.7 

T9* 

1.3 

1.2 

0.2 

1.0 

0.4 

T12 

T13 

0.2 

** 

- 

TU 

1.2 

0.15 

T15* 

0.6 

- 

T16 

0.1 

T17 

.1.4 

0.2 

t Total  range  of 

background  samples. 

0.2 

- 

NOTE: 


-  — x  W  associated  «tt  poUuted  a; 

nderlined  concentrations  are  greater  th* 

.  ,  greater  than  variation  in 

background  samples 


55 


I 


results  were  obtained  along  the  western  edge  of  the  sampling  grid 
(T4-T8),  but  positive  results  were  also  seen  along  the  centerline 
(T16  and  T13)*  These  observations  indicate  that  the  wind  direc¬ 
tion  probably  shifted  slightly  from  the  time  it  was  measured 
(13:00)  to  the  time  of  the  last  fuel  dump  (15:15)*  It  is 
interesting  that  some  stations  saw  concentrations  well  above  the 
expected  levels  while  adjacent  stations  saw  nothing.  This 
variation  is  probably  due  in  part  to  the  non-uniformity  of  the 
fuel  dump  plume  after  the  vortex  structure  begins  to  break-up  and 
fold  on  itself  under  the  influence  of  large  scale  atmospheric 
eddies.  The  overall  results  are  felt  to  be  in  general  agreement 
with  expectations,  demonstrating  that  only  minimal  concentrations 
of  fuel  vapor  are  reached  at  the  ground,  even  for  very  low  alti¬ 
tude  releases. 


56 


SECTION  IV 


PREDICTED  GROUNDFALL  OF  JETTISONED  JP-4 
Coupling  the  fuel  droplet  evaporation  and  free-fall  model  with  the 

h  : 

experimentally  determined  initial  size  distribution,  it  is  now 

possible  to  predict  the  fate  of  fuel  jettisoned  under  differing 

r 

circumstances.  In  all  cases  these  predictions  apply  only  for  Jet¬ 
tisoning  of  JP-4  fuel,  although  estimates  for  other  fuels  will  be 
provided  where  possible.  Atmospheric  dispersion,  which  Is  not 
considered  in  the  droplet  model,  will  also  be  treated  in  order  to 
fully  describe  the  physical  fate  of  the  fuel. 

4.1  FRACTION  OF  THE  FUEL  REACHING  THE  GROUND 

Figures  14  and  15  show  the  percentage  of  JP-4  fuel  which  can  be 
expected  to  reach  the  ground  before  evaporating,  based  on  the 
release  altitude  and  the  temperature  at  ground-level.  These  curves 
indicate  that  except  for  very  low  altitude  releases,  the  fraction 
of  fuel  reaching  the  ground  as  liquid  droplets  is  primarily  deter¬ 
mined  by  the  temperature.  For  JP-4  fuel  Jettisoned  higher  than 
1500  meters  and  at  temperatures  above  freezing  (0°C),  more  than  98 
percent  of  the  fuel  should  evaporate  before  reaching  the  ground. 

For  less  volatile  fuels,  such  as  the  Air  Force's  JP-8,  commercial 
Jet-A,  or  the  Navy’s  JP-5,  much  more  of  the  fuel  will  reach  the 
ground  as  liquid  droplets  than  would  be  the  case  for  JP-4  under 
the  same  conditions.  Preliminary  calculations  suggest  that  the 
percentage  of  JP-8,  Jet-A,  or  JP-5  reaching  the  ground  at,  for 
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example,  0°C  would  be  somewhat  between  the  curves  for  JP-4  at 
-20°C  and  -40°C  in  Figure  14. 

Figure  16  compares  the  fuel  droplet  size  distribution  for  the 
KC-135  ("Good  -  Corrected")  which  was  used  in  these  calculations, 
with  the  smallest  experiment  distribution  obtained  for  the 
Buccaneer  ("Cross  and  Picirnett  -  Experimental")  and  the  suggested 
larger  distribution  for  other  aircraft  ("Cross  and  Picknett  - 
Composite")  from  Reference  6.  These  distributions  probably 
represent  the  range  of  droplet  sizes  that  can  be  expected  for  most 
aircraft.  The  effect  of  these  different  distributions  on  the  pre¬ 
dicted  amount  of  liquid  fuel  reaching  the  ground  is  shown  in 
Figure  17.  Also  shown  are  the  predictions  for  a  single  droplet 
whose  diameter  is  equal  to  the  mass  median  diameter  of  the  KC-135 
distribijtion  (27h  microns).  The  difference  between  the  predic¬ 
tions  based  on  the  KC-135  distribution  and  those  obtained  with  the 
other  aircraft  distributions  is  generally  less  than  40  percent. 
That  is,  the  ground-fall  of  fuel  from  different  aircraft  should 
not  vary  more  than  0.6  to  1.4  times  the  estimates  shown  in  Figures 
14  and  15-  The  surprising  agreement  of  the  predictions  for  a 
single  droplet  at  the  mass  median  diameter  with  those  for  the 
whole  distribution  indicates  that  the  central  tendency  or  spread 
of  the  droplet  sizes  is  relatively  unimportant  for  determing  the 
composite  evaporation  and  free-fall  of  the  distribution. 

In  addition  to  experimentally  determining  an  initial  droplet  size 
distribution  for  fuel  jettisoning.  Cross  and  Picknett  also  esti¬ 
mated  the  fraction  of  fuel  reaching  the  ground  as  a  function  of 
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PERCENT  MASS  PER  MICRON 


Figure  16.  Experimental  Fuel  Droplet  Size  Distribution  for 
Aircraft  Fuel  Jettisoning 
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PREDICTED  LIQUID  FUEL  REACHING  THE 
OTHER  DISTRIBUTIONS  (PERCENT) 


PREDICTED  LIQUID  FUEL  REACHING  THE  GROUND  USING 
DISTRIBUTION  FOR  THE  KC-135  (PERCENT) 

Figure  17.  Effect  of  I.iitial  Droplet  Size  Distribution  on  the 
Percent  of  Fuel  Reaching  the  Ground 
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release  height  and  temperature  (Reference  6).  Their  estimates 
were  based  on  the  droplet  modelling  results  of  Lowell  (Reference 
2)  for  a  fall  of  300  meters.  Unfortunately,  to  model  release 
heights  greater  than  300  meters  they  calculated  the  evaporation 
during  successive  300  meter  intervals  by  assuming  that  the 
droplets  would  follow  the  same  curve  of  percent  evaporation  ver¬ 
sus  diameter  used  for  the  first  300  meters.  However,  as  discussed 
in  Section  2.2  of  this  report  the  evaporation  of  a  fuel  droplet 
changes  not  only  due  to  its  decreasing  diameter,  but  also  due  to 
its  changing  composition.  The  more  volatile  components  are 
stripped  away  in  the  early  evaporation,  leaving  only  the  slower 
evaporating  components  in  the  droplet.  By  ignoring  composition 
changes  Cross  and  Picknett  grossly  underestimated  the  fraction  of 
fuel  reaching  the  ground  for  release  heights  of  greater  than  300 
meters;  therefore,  their  predictions  cannot  be  used. 

4.2  OTHER  MODEL  PREDICTIONS 

As  can  be  seen  in  Figure  18,  the  vast  majority  of  the  fuel  evap¬ 
orates  in  the  first  few  minutes  after  release.  For  fuel  released 
at  1500  meters  when  the  ground-level  temperature  is  above  freezing, 
less  than  10  percent  of  the  fuel  remains  after  10  minutes.  During 
this  short  time  interval  even  the  largest  droplet'0  only  fall 
700-900  meters,  and  most  of  the  vapor  is  generated  within  a  few 
hundred  meters  of  the  release  height.  The  implication  of  these 
calculations  for  dispersion  modeling  is  that  the  fuel  dump  can  be 
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Figure  18.  Evaporation  of  a  Distribution  of  Droplets 
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approximated  as  an  elevated  vapor  source  at  the  release  height 
without  seriously  degrading  the  accuracy  of  the  estimated  ground- 
level  concentrations.  This  approximation  was  used  for  the  disper¬ 
sion  calculations  in  this  report. 

Due  to  the  different  settling  velocities  of  droplets  with  varying 
diameters,  not  all  of  the  liquid  fuel  reaches  the  ground  at  the 
same  time.  Two  examples  of  the  time  spread  over  which  liquid  fuel 
reaches  the  ground  from  1500  meters  are  shown  in  Figure  19.  At  a 
ground-level  temperature  of  -20°C,  10  percent  of  the  fuel  reaches 
the  ground  in  liquid  form.  The  largest  droplets,  on  the  order  of 
300  microns,  reach  the  ground  in  less  than  20  minutes.  Most  of 
the  unvaporized  fuel  reaches  the  ground  during  the  first  hour 
after  release,  but  a  significant  fraction  takes  much  longer. 

Raising  the  ground-level  temperature  reduces  the  fraction  of  fuel 
reaching  the  ground  in  liquid  form,  and  also  delays  the  onset  of 
ground-fall.  Although  increasing  the  release  altitude  above  1500 
meters  does  not  substantially  lower  the  fraction  of  liquid  fuel 
reaching  the  ground,  it  greatly  increases  the  time  required  for 
the  droplets  to  reach  the  ground.  The  time  of  fall  for  the  first 
droplets  to  reach  the  ground  at  several  altitudes  and  temperatures 
is  shown  in  Table  12.  These  predictions  of  the  fuel  droplet  model 
can  be  approximated  by  the  formula  T  =  NH,  where  T  is  the  time  of 
fall  in  minutes,  H  is  the  release  height  in  kilometers,  and  N  is 
18,  12,  and  10  minutes  per  kilometer  for  ground-level  temperatures 
of  0°C,  -20°C,  and  -40°C,  respectively.  These  results  can  be  used 
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to  estimate  the  effect  of  winnowing  by  the  wind  on  the  dispersion 
of  the  fuel. 

For  a  ground-level  temperature  of  -20°C,  slightly  over  10  percent 
of  the  fuel  released  at  1500  meters  will  reach  the  ground  as 
droplets.  The  different  size  droplets  will  be  distributed  in 
time,  as  indicated  by  the  upper  scale  in  Figure  19»  This  same 
distribution  is  shown  on  a  linear  scale  in  Figure  20.  The  modal 
diameter  at  -20°C  is  190  microns.  For  a  ground-level  temperature 
of  0°C  evaporation  is  more  complete,  and  the  modal  diameter  at 
the  ground  is  less  than  100  microns. 

The  composition  of  the  fuel  droplets  which  reach  the  ground  is  no 
longer  the  same  as  that  of  the  JP-4  fuel  which  was  jettisoned. 

The  more  volatile,  lower  molecular  weight  components  will  have 
evaporated  off,  leaving  a  residual  mixture  of  the  higher  molecular 
weight  components.  A  typical  composition  for  a  distribution  of 
droplets  which  has  evaporated  to  10  percent  of  the  original  mass 
is  shown  in  Table  13*  As  can  be  seen  in  this  table,  the  final 
droplet  composition  resembles  kerosene  rather  than  JP-4,  in  that 
it  is  composed  chiefly  of  hydrocarbons  having  a  molecular  weight 
greater  than  150  grams  per  mole. 

4.3  DISPERSION  OF  THE  FUEL 

Once  the  fraction  of  fuel  reaching  the  ground  in  liquid  form  has 
been  determined,  the  fate  of  the  liquid  and  vapor  fractions  can 
be  Investigated  separately.  The  next  step  in  assessing  the  envir¬ 
onmental  consequences  of  a  jettisoning  incident  is  to  estimate  the 


TABLE  13.  RESIDUAL  COMPOSITION  OP  FUEL  DROPLETS 
AFTER  90  PERCENT  OF  THE  ORIGINAL  MASS  HAS  EVAPORATED 


Original  Percent  Percent  Percent 

Percent  of  Remaining  of  of  Initial  of  Final 
Droplet  Mass  Component  Droplet  Mass  Droplet  Mass 


iso-pentane 

3.2 

0 

- 

- 

iso-hexane 

7.1 

0 

- 

- 

cyclohexane 

2.2 

0 

- 

- 

benzene 

0.3 

0 

- 

- 

3-methylhexane 

8.6 

0 

- 

- 

methyl cyclohexane 

7.3 

0 

- 

- 

toluene 

0.8 

0 

- 

- 

4-methylheptane 

9.4 

0 

- 

- 

cis-1,  4-dimethyl- 
cyclohexane 

_ 7 . 7 

0 

. 

— . —  '  - _ 

m-xylene 

1.8 

0 

- 

4-methylootane 

8.7 

0 

- 

- 

isopropyl- 

cyclohexane 

4.6 

0.3 

0.01 

0.1 

l-ethyl-2  methyl- 
benzene 

2.8 

1.6 

0.05 

0.5 

2, 7-dimethyloctane 

7.0 

1.0 

0.07 

0.7 

p-methane  (cis) 

3.9 

4.2 

0.17 

1.6 

p-cymene 

2.1 

7.4 

0.15 

1.5 

napthalene 

0.3 

46.4 

0.13 

1.3 

undecane 

4.7 

26.9 

1.29 

12.6 

3-methylbutyl- 

cyclohexane 

2.7 

27.0 

0.73 

7.1 

3-methylenedecalin 

(trans) 

4.0 

31.0 

1.26 

12.3 

1-bu ty 1-3-methyl - 
benzene 

1.2 

35.4 

0.45 

4.4 

1-methylnapthalene 

0.3 

66.8 

0.18 

1.7 

dodecane 

2.8 

48.3 

1.36 

13.3 

3-ethylbutyl- 

cyclohexane 

1.3 

42.9 

0.55 

5.4 

1,  3,5-triethyl¬ 
benzene 

0.6 

47.0 

0.27 

2.6 

2, 3-dimethyl- 
napthalene 

0.3 

78.9 

0.21 

2.0 

tridecane 

1.1 

63.2 

0.71 

6.9 

3-isopropylbutyl- 

cyclohexane 

0.4 

55.7 

0.24 

2.3 

3, 5-diethyl-l- 
propylbenzene 

0.1 

61.6 

0.07 

0.7 

tetradecane 

0.2 

73.8 

0.15 

1.5 

pentadecane 

0.1 

81.4 

0.08 

C.8 

perhydro- 

phenanthrene 

2.2 

86.7 

1.96 

19.1 

residual 

0.2 

99.8 

0.16 

1.6 
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extent  to  which  the  fuel  vapor  and  droplets  are  dispersed  prior  to 
reaching  ground-level. 

4.3.1  VAPOR  DISPERSION 

Since  most  of  the  fuel  jettisoned  from  an  aircraft  will  evaporate 
within  a  few  hundred  meters  of  the  release  height,  the  concentra¬ 
tion  of  hydrocarbons  at  ground-level  is  primarily  determined  by 

the  ability  of  atmospheric  diffusion  processes  to  transport  the 

,  / 

vapors  downward.  The  highest  ground-level  concentrations  will  be 
obtained  under  unstable  atmospheric  conditions,  since  the  vertical 
transport  is  then  a  maximum.  The  longest  exposures  will  occur 
when  the  fuel  is  jettisoned  parallel  to  the  wind. 

The  gradient-transfer  (K)  theory  of  atmospheric  diffusion  can  be 

used  to  estimate  the  "worst-case”  vapor  concentrations  at  the 

ground.  According  to  K  theory,  the  vertical  spread  of  a  plume 

(as  measured  by  its  standard  deviation,  0  )  can  be  related  to  the 

z 

vertical  component  of  the  atmospheric  eddy  diffusivity;  K  ,  by  the 

z 

2 

formula  K„  =  /2T,  where  T  is  the  time  since  release  (Reference 

11)  .  A  profile  of  eddy  diffusivity  in  the  mixing  layer  (Reference 

12)  shows  that  K_  is  generally  less  than  100  square  meters  per 
second,  even  under  relatively  unstable  conditions.  For  a  Gaussian 
plume  profile,  the  maximum  ground  concentration  Is  reached  when 

a„  =  0.70TH,  where  H  is  the  release  height.  Using  a  value  for  K 
of  100  square  meters  per  second  to  represent  strong  vertical  mix¬ 
ing,  the  time  at  which  the  maximum  ground  concentration  is  reached 
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2  2 

is  given  by  Tmax  =  o  /2K  =  H  /40Q.  Expressing  T  in  minutes  and 

Ci  Cl 

p 

H  in  kilometers,  this  formula  becomes  T  =  40H  . 

The  horizontal  spread  (c  )  of  the  plume  during  this  period  can  be 

J 

estimated  by  o  =  0.5t  (Reference  13)*  This  value  provides  a 
«y 

good  fit  to  experimental  data  on  a  for  plume  travel  periods  from 

V 

1-100  hours  (Reference  11).  Assuming  that  the  highest  con¬ 
centrations  in  the  plume  occur  over  a  width  V/  =  2cv,  the  width  of 
the  plume  at  the  time  the  maximum  ground-level  concentrations  are 
reached  is  W  =  H^/400.  For  W  and  H  in  kilometers,  W  =  2.4H^. 

Using  a  simple  box  model  and  ignoring  upward  diffusion  of  the 
plume,  we  can  then  conservatively  estimate  the  ground  concentration 
with  the  formula  X  =  1000Q/VWH,  where  X  is  the  concentration  in 
micrograms  per  cubic  meter,  Q  is  the  jettison  rate  in  kilograms 
per  second,  and  V  is  the  aircraft  velocity  in  meters  per  second. 
Substituting  for  W  from  the  formula  derived  above  yields 
X  =  400Q/VH  .  The  predications  of  this  box-model  for  several  fuel 
jettisoning  scenarios  are  shown  in  Table  14. 

TABLE  14.  WORST-CASE  GROUND-LEVEL  VAPOR  CONCENTRATIONS 


Release 

Aircraft 

Jettison 

Maximum  Ground 

-Level 

Height 

Velocity 

Rate 

Concentration 

(ug/nw) 

(kilometers) 

(m/s) 

(kg/s) 

K-Theory/Box  Gaussian 

.3 

175 

17 

1439 

298 

.775 

175 

56 

275 

147 

1.5 

175 

5 

3.4 

3.5 

6 

175 

56 

0.6 

0.4 

An  alternative  technique  for  estimating  the  worst-case  con¬ 
centrations  of  vapor  at  ground-level  is  to  rely  solely  on  a 
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Gaussian  diffusion  analysis  (Reference  10).  In  this  case  the 
maximum  ground-level  concentration  is  given  (as  in  Section  3*4.1) 
by  X  -  Q/fTOyCTzUe,  where  normally  U  is  the  windspeed,  but  here 
U  =  V.  Then  assuming  a.,  =  c„  (true  within  a  factor  of  two  for 

J  6 

2 

moderately  unstable  atmospheric  conditions) ,  X  =  2Q/7reVH  .  For  X 
in  micrograms  per  cubic  meter,  Q  in  kilograms  per  second,  and  H 

p 

in  kilometers,  this  formula  bedcomes  X  =  276Q/VH  .  The  predic¬ 
tions  of  this  Gaussian  model  are  also  shown  in  Table  14.  The  two 
models  agree  for  a  release  height  of  1.5  kilometers  and  diverge 
elsewhere.  The  theoretical  basis  for  the  two  models  is  quite 
different;  the  Gaussian  model  is  most  valid  for  low  altitudes, 
whereas  the  K-theory  box-model  applies  best  to  high  altitudes 
near  the  top  of  the  mixing  layer.  Thus  the  two  analyses  are 
complimentary.  As  a  means  of  comparison  with  a  more  exacting 
model,  the  Environmental  Protection  Agency's  "Point-Area-Line" 
(PAL)  computer  dispersion  model  (Reference  14)  was  run  for  the 
case  in  Table  14  for  1.5  kilometers.  Under  moderately  unstable 
conditions  (stability  category  B)  with  a  windspeed  of  2  meters 
per  second,  the  highest  ground-level  concentration  predicted  by 
the  PAL  model  was  8.5  micrograms  per  cubic  meter,  in  reasonable 
agreement  with  the  two  simplified  approximation  techniques. 


4.3.2  DROPLET  DISPERSION 


At  temperatures  below  freezing  the  fraction  of  fuel  reaching  the 
ground  as  liquid  droplets  can  become  significant.  As  in  the  case 
of  fuel  vapor,  however,  natural  atmospheric  dispersion  would 
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reduce  the  droplet  density  at  the  ground.  For  jettisoning 
parallel  to  the  wind,  the  droplets  would  be  spread  over  nearly  as 

p 

wide  a  path  as  the  vapor:  W  =  2.4HC,  where  W  and  H  are  in 
kilometers.  Therefore,  the  maximum  liquid  fuel  contamination  of 
the  ground,  C,  (in  milligrams  per  square  meter)  would  be  given 
approximately  by  the  formula  C  =  10PQ/(VW),  where  P  is  the  per¬ 
centage  of  fuel  reaching  the  ground  in  liquid  form  from  Figure  14 
or  15. 

For  jettisoning  in  any  direction  other  than  parallel  to  the  wind, 
the  droplets  would  be  further  separated  by  the  process  of  win¬ 
nowing:  as  the  dr>  /le4‘  are  carried  along  by  the  wind  the  larger 
droplets  fall  fs •  ?  nd  are  deposited  sooner  (closer  to  the  jet¬ 
tisoning  locvt. ’  .;n;  than  the  smaller  droplets,  which  tend  to  be 

.;v. 

transported  more  like  the  vapor.  The  distance  in  kilometers  that 

the  fuel  will  be  carried  downwind  before  reaching  the  ground  is 

given  by  L  =  0.06UT,  where  U  is  the  average  windspeed  (in  meters 
per  second)  between  the  jettisoning  altitude  and  the  surface,  and 

p 

T  is  the  time  of  fall  in  minutes.  As  stated  earlier,  T  =  40H  for 

the  vapor,  while  for  the  largest  droplets  T  =  NH,  where  N  is  18, 

12  and  10  minutes  per  kilometer  at  0°C,  -20°C,  and  -40°C, 
respectively.  Therefore,  for  fuel  jettisoned  crosswind,  the 
droplets  will  be  spread  by  winnowing  over  a  distance 

p 

W  =  0.06U(40H  -  NH) .  Several  cases  comparing  jettisoning 

parallel  and  perpendicular  to  the  wind  are  shown  in  Table  15. 
Generally,  the  average  liquid-fuel  ground  contamination  is 
reduced  by  a  factor  of  three  to  four  when  jettisoning  crosswind. 
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“"  *-"  ^AAMELES  OP  LIQUID  PIIRT 
CONTAMINATION  OF  THE  GROUND 


Aircraft : 

Airspeed  (m/s): 
Jettison  Rate  (kg/s): 
Release  Height  (km) : 
Ground-Level 
Temperature  (c) : 
Windspeed  (m/s) : 


Parallel  to  Wind: 
Width  at  Ground  (km): 
Concentration  (mg/m 2)  • 


Perpendicular  to  Wind: 
Width  at  Ground  (km): 
Concentration  (mg/m 2) ■ 


F-i| 

P-111 

FB-I11 

KC-135 

P-111 

175 

5 

1.5 

175 

17 

1.5 

175 

17 

6 

175 

50 

6 

175 

17 

.3 

0 

5 

0 

5 

-20 

3 

-20 

11 

20 

A 

5 

.08 

5 

.28 

90 

.09 

90 

.27 

.2 

3.1 

19 

.02 

19 

.08 

250 
•  03 

330 

.07 

.2 

3.1 
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SECTION  V 


SUMMARY  AND  CONCLUSIONS 

A  computer  model  has  been  developed  which  simulates  the  evapora-, 
tion  and  free-fall  of  fuel  droplets  in  the  atmosphere.  This  model 
gives  good  agreement  with  evaporation  rates  of  fuel  droplets 
measured  experimentally.  A  study  of  fuel  jettisoning  by  a  KC-135 
tanker  aircraft  was  performed  to  obtain  a  fuel  droplet  size  dis¬ 
tribution  for  input  to  the  model.  The  droplet  size  distribution 
was  measured  directly  in  flight,  but  only  after  significant  evap¬ 
oration  was  known  to  have  occurred.  The  observed  droplet  sizes 
were  therefore  corrected  for  evaporation  using  the  fuel  droplet 
evaporation  model.  Correcting  the  observed  droplet  mass  densities 
for  evaporation,  all  of  the  jettisoned  fuel  was  accounted  for. 

The  original  distribution  thus  obtained  has  a  mass  median  diameter 
of  265  microns  with  a  maximum  droplet  size  of  around  500  microns. 
During  this  study,  sampling  was  also  performed  at  ground  level  to 
determine  whether  the  jettisoned  fuel  reached  the  ground  in  signi¬ 
ficant  concentrations.  For  fuel  jettisoned  as  low  as  750  meters 
above  the  ground  at  temperatures  around  11°C,  no  liquid  fuel  could 
be  detected  by  ground  observers  and  no  significant  hydrocarbon 
concentrations  (greater  than  a  few  ppmC)  were  measured  by  the 
sampling. 

Based  on  the  fuel  droplet  evaporation  and  free-fall  model,  more 
than  98  percent  of  JP-^t  fuel  will  evaporate  before  reaching  the 
ground  if  jettisoned  higher  than  1500  meters  when  the  ground-level 
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temperature  is  above  freezing.  The  droplets  and  vapor  will  be 
widely  dispersed  by  atmospheric  turbulence,  quickly  resulting  in 
a  fuel  density  too  low  to  create  any  perceptible  environmental 
changes.  The  formulas  and  graphs  in  Section  IV  can  be  used  to 
estimate  the  likelihood  of  significant  ground-level  concentrations 
of  fuel  vapor  or  liquid  following  a  specific  fuel  jettisoning 
incident. 

This  report  has  dealt  only  with  the  short-term  physical  fate  of 
the  jettisoned  fuel.  The  possible  longer-term  fates  of  the  fuel 
vapor  in  the  atmosphere  are  discussed  in  Reference  1.  Since  the 
vapor  is  quickly  dispersed  and  diluted  below  the  levels  at  which 
it  could  be  harmful  in  itself  (around  500  ppm),  its  environmental 
impact  derives  principally  from  its  role  in  the  production  of 
photochemical  oxidant  pollution  (ozone  and  smog).  Studies  of  this 
role  are  currently  in  progress.  The  possible  long-term  fate  of 
any  liquid  fuel  which  reaches  the  ground  is  also  under  study,  with 
particular  emphasis  on  the  chemical  and  biological  fate  of  the 
hydrocarbons  in  an  aqueous  environment.  When  these  studies  are 
completed  the  overall  environmental  consequences  of  fuel  jet¬ 
tisoning  can  be  determined. 
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APPENDIX  A 

EFFECT  OF  DROPLET  TEMPERATURE  CALCULATIONS 
IN  THE  FUEL  DROPLET  EVAPORATION  MODEL 

The  evaporation  rate  of  a  fuel  droplet  depends  in  part  on  the 
vapor  pressure  of  its  components  at  the  droplet  temperature.  As 
described  in  Appendix  D,  the  fuel  droplet  model  performs  an  energy 
balance  to  derive  the  steady  state  droplet  temperature.  This 
energy  balance  considers  heat  transfer  by  conduction  with  the  air, 
absorb tion  and  emission  of  radiation,  solar  insolation,  and  evap¬ 
orative  cooling.  In  addition  the  model  calculates  an  initial 
droplet  temperature  higher  than  that  of  the  ambient  air,  assuming 
equilibration  of  the  fuel  in  the  tank  with  the  skin  of  the  air¬ 
craft  before  release.  The  model  could  be  greatly  simplified  by 
eliminating  these  calculations  and  assuming  instead  that  the 
droplet  temperature  is  always  equal  to  to  the  temperature  of  the 
ambient  air.  Therefore,  it  is  worthwhile  to  investigate  the 
importance  of  the  energy  balance  and  initial  droplet  temperature 
routines . 

The  effect  of  these  two  calculations  on  the  droplet  temperature  is 
shown  in  Figure  A-l.  Due  to  the  dominant  effect  of  evaporative 
cooling  in  the  energy  balance,  the  droplet  cools  appreciably  below 
the  ambient  air  temperature  uuring  the  initial  moments  when  the 
highly  volatile  components  are  being  stripped  away.  However,  the 
combined  effects  of  solar  insolation  and  heat  conduction  from  the 
air  rapidly  restore  the  droplet  to  the  ambient  temperature  when 
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GROUND- LEVEL  TEMPERATURE;  t-20°C 
INITIAL  ALTITUDE  :  1500  METERS 
INITIAL  DROPLET  DIAMETER  :  270  MICRONS 


-WITH  ENERGY  BALANCE  AND  INITIAL  DROPLET 
TEMPERATURE  ABOVE  AMBIENT 


NO  ENERGY  BALANCE 


balance 


WITH  ENERGY  BALANCE 


.01  .1 
TIME  (MINUTES) 


Figure  A-l.  Difference  Between  Droplet  and  Ambient  Temperatures 
for  Different  Calculations 


GROUND- LEVEL  TEMPERATURE:  +  20°C 
INITIAL  ALTITUDE;  1500  METERS 
INITIAL  DROPLET  DIAMETER  :  270  MICRONS 


•WITH  ENERGY  BALANCE  AND  INITIAL  DROPLET 
TEMPERATURE  ABOVE  AMBIENT 


•  NO  ENERGY  BALANCE 


\  / —  WITH  ENERGY  BALANCE 


.01  .1 
TIME  (MINUTES) 


Figure  A-2.  Predicted  Evaporation  Rates  for  Different  Droplet 
Temperature  Calculations 


evaporation  slows.  The  entire  droplet  cooling  and  warming  cycle 
lasts  less  than  a  minute  in  most  cases.  The  assumption  of  an  ini¬ 
tial  droplet  temperature  above  ambient  tends  to  delay  and  partially 
offset  the  maximum  cooling  effect,  but  does  not  eliminate  it. 

Figure  A-2  shows  the  resulting  variation  in  the  droplet  evapora¬ 
tion  rates.  Although  assuming  an  initial  droplet  temperature 
above  ambient  produces  a  larger  change  in  instantaneous  evapora¬ 
tion  rate,  the  opposing  effect  of  the  energy  balance  is  of  longer 
duration.  The  net  effect  on  the  predicted  mass  remaining  in  the 
droplet  as  a  function  of  time  is  shown  in  Figure  A-3.  The  raised 
initial  temperature  causes  the  droplet  mass  to  fall  off  much  more 
rapidly  at  first,  but  evaporative  cooling  quickly  counters  this 
effect,  and  after  the  first  few  seconds  the  droplet  mass  is 
greater  than  it  would  be  if  both  effects  were  ignored.  Regard¬ 
less  of  the  calculations  used,  the  droplet  mass  remaining  is 
essentially  the  same  after  the  first  minute.  The  reason  the 
three  different  calculations  come  back  together  after  initially 
diverging  is  that  the  changing  composition  of  the  droplet  is  a 
factor  in  the  evaporation  rate.  For  example  if  the  droplet  is 
made  to  evaporate  more  slowly,  its  composition  will  include  more 
of  the  volatile  components.  Then  when  the  restraint  on  evapora¬ 
tion  is  removed  the  more  volatile  composition  of  the  droplet  will 
accelerate  its  evaporation,  counteracting  the  effect  of  the 
earlier  restraint. 

The  primary  conclusion  which  can  be  drawn  from  these  results 
is  that  for  predictions  of  droplet  evaporation  over  periods  greater 


82 


than  one  minute,  the  droplet  temperature  can  be  assumed  to  be 
identical  to  the  ambient  air  temperature.  Only  for  predictions  of 
the  droplet  evaporation  rate  and  history  during  the  first  minute 
are  the  energy  balance  and  initial  droplet  temperature  calculation 
important. 
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APPENDIX  R 


DETAILED  DESCRIPTION  OP  THE  FUEL  DROPLET  MODEL 

The  fuel  droplet  evaporation  and  free-fall  model  breaks  a  given 
droplet’s  fall  into  a  series  of  small  time  intervals.  The  dis¬ 
tance  of  fall  during  each  interval  is  calculated  assuming  the 
droplet  is  falling  at  the  terminal  velocity  for  its  current 
diameter,  density  and  altitude.  Loss  of  mass  through  evaporation 
is  calculated  assuming  Raoult’s  law;  that  is,  each  component  evap¬ 
orates  independently.  An  energy  balance  routine  adjusts  the 
droplet  temperature  to  allow  for  evaporative  cooling,  radiation, 
conduction  and  insolation  effects.  The  new  droplet  composition, 
mass,  and  altitude  are  used  as  initial  conditions  for  the  next 
interval.  This  stepwise  approximation  continues  until  the  droplet 
impacts  on  the  ground  or  loses  99*9  percent  of  its  Initial  mass. 

The  Initial  conditions  which  must  be  known  are  the  droplet's  ori¬ 
ginal  composition,  altitude  and  diameter;  the  temperature  and 
altitude  at  local  ground  level;  and  the  aircraft's  air  speed.  The 
initial  droplet  temperature  is  then  taken  as  the  corresponding 
stagnation  temperature,  assuming  equilibration  of  the  fuel  tanks 
with  the  skin  of  the  aircraft.  In  the  early  intervals  the  droplet 
is  allowed  to  cool  through  evaporation  until  an  energy  balance  is 
achieved. 

At  the  beginning  of  each  interval,  the  droplet's  current  mass  com¬ 
position  and  altitude  are  known.  (For  the  first  interval,  the  mass 
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is  determined  from  the  droplet’s  diameter  and  density;  calculation 
of  the  latter  is  discussed  below).  The  lo<;al  air  temperature  is 
derived  from  the  droplet  altitude  and  the  temperatur*  v  at' the 
ground  using  the  standard  lapse  rate  of  6.5  K/km.  Atmospheric 
pressure,  density,  and  viscosity  are  then  calculated  as  recommended 
in  the  US  Standard  Atmosphere,  1976  (Reference  B-l) : 

P  =  P  (T  /T  )  (£Ma/R0L) 
a  o  '  o'  a' 


where  P.  =  the  local  atmospheric  pressure  (N/m') 

cL 

2 

P0  =  atmospheric  pressure  at  sea  level  =  101325  N/m 
Tq  =  sea  level  temperature  (K) 

Tq  =  local  air  temperature  (K) 

a 

2 

g  =  acceleration  of  gravity  =  9.81  m/s 
Ma  =  molecular  weight  of  air  =  28.96  kg/kmol 
Rq  =  universal  gas  constant  =  831^  N*m/K‘kmol 
L  =  standard  lapse  rate  =  0.0065  K/m 


Po  =  P0M0/.RrtT<k 
Ka  a  a  o  a 


3  v 

where  p  =  the  local  air  density  (kg/mJ) 


1.458  x  10“6  Tq  1,5 

a 

U  =  - 

a  110.4  +  T 

a 

2 

where  =  the  local  air  viscosity  (N.s/m  ) 

For  the  first  interval  the  droplet  temperature  is  assumed  to  be 
at  the  equilibrium  temperature  of  the  aircraft  fuel  tanks,  which 
is  approximately  the  stagnation  temperature  corresponding  to  the 
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aircraft  velocity  and  local  air  temperature  (Reference  B-2).  Tile 
stagnation  temperature,  Tg  (K) ,  can  be  calculated  from  the  -true 
air  speed,  V  (m/s),  and  the  local  air  temperature:,,  .  ,  . 

Ts  -  Ta  (1  +  V2/5Cs2) 

where  CH  is  the  local  speed  of  sound  (m/s)  calculated  as  in 
Reference  B-l: 

°s  ”  tl-*»  R0VKa)V2 

o 

The  density,  (kg/m  )  is  calculated  for  each  component  in  the 
mixture  using  a  linear  expansion  coefficient: 

P±  =  20Pi/(l  +  0.001  (Td  -  293.15)) 

20  3 

where  p1  =  the  component’s  density  at  20C  (kg/irr) 

Td  =  the  droplet  temperature  (K) 

An  expansion  coefficient  of  0.001  was  chosen  to  give  the  best  fit 
for  a  variety  of  hydrocarbons.  Assuming  an  ideal  solution,  the 
droplet’s  volume  is  then  a  sum  of  its  components. 


The  droplet  terminal  velocity,  Ut(m/s)  is  determined  following 
the  method  of  Lowell  (References  B-3  and  B-4) :  A  parameter  "q” 

1/p 

(equivalent  to  Lowell’s  <}>  c)  is  defined: 

q  =  ReC^2  =  (4PaPdgd3/3na?)  ^2 

where  Re  =  the  Reynolds  number  of  the  system 
CD  =  the  drag  coefficient 
p,  =  the  droplet  density  (kg/m3) 


d  =  the  droplet  diameter  (m) 

An  empirical  relation  between  q  and  Re  is  then  used  to  find  Re: 

In  Re  =  -3.13  +  2.06  In  q  -  0.083(ln  q)2 

This  relation  was  chosen  to  fit  the  AEDC  data  (Reference  B-5)  for 
Reynolds  numbers  over  100,  while  agreeing  with  Stoke' s  law 
(Re  =  q2/2iO  for  Re  <  1.  A  plot  showing  the  success  of  this  fit 
is  presented  in  Figure  B-l.  The  droplet  terminal  velocity  is 
then : 


Ut  =  >‘aRe/d‘,a 

and  the  distance  the  droplet  falls  during  the  interval  is 
approximated  by  the  product  of  Ut  with  the  duration  of  the 
interval,  At( s) , 

The  evaporation  of  a  falling  single-component  drop  is  treated  in 
Bird,  Stewart,  and  Lightfoot  (Reference  B-6).  To  treat  the 
multicomponent  fuel  droplet  the  contribution  from  each  component 
is  calculated  separately  assuming  independent  transport  and 
ignoring  any  internal  resistance: 

2 

Am^  =  ird  hj_Pj_E  jA  t 

where  Am.  =  the  mass  of  component  "i"  evaporated  during  the 
“  time  At  (kg) 

h,  =  the  mass  transfer  coefficient  for  the  component 
( s/m) 

2 

p^  =  the  true  vapor  pressure  of  the  component  (N/m  ) 

Ej  =  the  mole  fraction  of  the  component  in  the 
x  droplet 


Figure  B-l.  Relation  Between  the  Reynolds  Number  and  the  Drag  Coefficient 


The  component  vapor  pressure  is  calculated  using  a  modified 
Antoine  equation: 


In  p±  =  20.53  -  2899/(385.15  VTb,i  ’  62,3) 
where  Tb  ^  =  the  normal  boiling  point  of  component  "i"  (K) 
The  constants  is  this  equation  were  obtained  by  using 
3, 3-dimethylhexane  as  a  reference  compound.  This  single  Antoine 
equation  can  be  used  for  all  of  the  components  through  the  use  of 
the  Ramsey-Young  Relation,  which  states  that  for  similar  substan¬ 
ces  (e.g.  hydrocarbons)  the  ratio  of  the  temperatures  at  which 
the  substances  exert  equal  vapor  pressures  is  the  same  for  any 
vapor  pressure.  As  used  in  the  Antoine  equation,  the  relation 
takes  the  form: 

Tr  "  Tb,r(Td/Tb,l) 

where  Tp  is  the  temperature  at  which  the  reference  compound 
(3»3-dlmethylhexane)  will  have  the  same  vapor  pressure  that  com¬ 
ponent  "i"  has  at  the  droplet  temperature,  and  p  and  Tb  1  are 
the  normal  boiling  points  of  the  reference  compound  and  component 
"i",  respectively.  Table  B-l  compares  the  vapor  pressure  pre¬ 
dicted  in  this  manner  with  experimental  vapor  pressures  for 
several  fuel  components  (from  Reference  B-7).  The  agreement  is 
generally  within  10  percent. 

The  mass  coefficient  is  a  function  of  the  component's  diffusivity 
2 

D  (m  /a),  the  Schmidt  number,  Sc,  and  the  Nusselt  number,  Nu: 


TABLE  B-l.  COMPARISON  OF  EXPERMENTAL 
AND  CALCULATED  VAPOR  PRESSURES 


Compound 


Vapor  Pressure  (Pa) 

Temperature (K)  Experimental  Calculated  Calc/Exp 


n-pentane 

223 

1333 

1197 

0.90 

n-pentane 

261 

13332 

13468 

1.01 

n-pentane 

309 

101325 

103839 

1.02 

iso-hexane 

214 

133 

109 

0.82 

iso-hexane 

243 

1333 

124  L 

0.93 

iso-hexane 

284 

13332 

13594 

1.02 

cyclohexane 

228 

133 

138 

1.04 

cyclohexane 

257 

1333 

1326 

0.99 

cyclohexane 

299 

13332 

13462 

1.01 

benzene 

236 

133 

292 

2.20 

benzene 

262 

1333 

1925 

1.44 

benzene 

299 

13332 

13990 

1.05 

2-methylhexane 

233 

133 

128 

0.96 

2-methylhexane 

264 

1333 

1367 

1.03 

2-methylhexane 

307 

13332 

13702 

1.03 

methyl cyclohexane 

237 

133 

98 

0.74 

methylcyclohexane 

270 

1333 

1204 

0.90 

toluene 

246 

133 

123 

0.92 

toluene 

280 

1333 

1429 

1.07 

4-methylheptane 

253 

133 

151 

1.14 

^ -methyl heptane 

285 

1333 

1420 

1.07 

dimethylcyclohexane 

254 

133 

113 

0.85 

dimethylcyclohexane 

288 

1333 

1253 

0.94 

m-xylene 

266 

133 

145 

1.09 

2, 7-dimethyloctane 

279 

133 

139 

1.05 

cis-decalin 

296 

133 

94 

diethylbenzene 

294 

133 

147 

0.71 

napthaline 

326 

133 

255 

1.92 

undecane 

306 

133 

180 

1.35 

triethylbenzene 

321 

133 

184 

1.38 

tridecane 

333 

133 

206 

1.55 

isopropylnapthal^ne 

349 

133 

154 

1.16 
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>5 


2.66  x  10“ JT 


-5m  2 


D1  * 


1  +  1 


M 


1 


M 


a 


Pa  (Vb1/3  +  0. 31) 2 


where  Tp  =  the  temperature  of  the  evaporating  film  (K) 

«  l/2(Ta  +  Td) 

=  the  molecular  weight  of  compound  "1"  (kg/kmol) 

=  the  molar  volume  of  compound  "1”  at  Its  normal 
boiling  point 

-  (1  +  0.001  (Tb>i  -  293.15)  )/20P1 


Nuj_ 

^1 


2  +  0.6  Re1/2  Sci1/3 
NulniMl/(dRoTa) 


The  equation  for  the  diffuslvlty  Is  derived  from  the  one  given  in 

Bird,  Stewart,  and  Lightfoot  by  using  a  typical  value  of  the 

Lennard-Jones  function,  0,  for  hydrocarbon-air  pairs  at  273K  and 

-1/2 

including  the  functional  dependence  of  n  on  Tp  (in  the  tem¬ 
perature  range  of  interest  )  explicitly.  The  error  introduced  by 
these  approximations  is  estimated  to  be  less  than  10  percent.  A 
comparison  of  experimental  and  predicted  dif fusivities  for 
several  pure  hydrocarbons  is  shown  in  Table  B-2. 

As  the  droplet  falls  it  cools  through  evaporation  and  heat 
transfer  until  it  is  accommodated  to  the  local  air  temperature. 
The  mass  lost  through  evaporation  is  included  in  a  total  energy 
balance  equation  (Reference  B-8)  to  determine  the  new  droplet 
temperature : 
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TflabAND‘MLoSMPmSMFpSSIE^PSNTAL 


Compound 


n-hexane 
2, 3-dime thylbutane 
methyl cyclopentane 
cyclohexane 
cyclohexane 
benzene 
toluene 
toluene 
n-octane 
n-octane 

f  th  *5“tr,imethylpentane 

ethylbenzene 

mesitylene 

propylbenzene 

n-decane 

napthalene 

dodecane 

diphenyl 


288 

288 

7.6 

7  C 

6.9 

288 

( •  5 

7  4. 

7.0 

288 

r .  o 

7 

7.4 

318 

1  •  0 

A  A 

7.5 

273 

0.0 

7  7 

9.0 

273 

303 

7.6 

8.8 

A  n 

7.2 

6.5 

27  3 

8.0 

303 

303 

273 

7.1 

7.1 

6  fi 

5.3 

6.6 

6.6 

273 

5.6 

6.0 

273 

5.5 

363 

h*  y 

8.4 

c  1 

5.6 

273 

8.4 

399 

J  •  A 

3.1 

6.1 

5.9 

273 

9.2 

5.0 

0.91 

0.93 

0.97 

0.99 

1.05 

0.94 

0.86 

0.91 

1.06 

0.93 

0.93 

0.91 

0.98 

1.14 

1.00 

1.16 

1.14 

0.82 


Q  *  TTd2At(eaaTa4  +  i/4(l-a)L  +  k(Ta-Td)  -  c^oT^)  -  Ml^m 
where  e&  *  the  emisslvlty  of  air  *  0.75 

ed  *  the  emisslvlty  of  tie  droplet  *  0.95 
o  *  the  Stef an-Boltzmann  constant 
*  5.67  x  10"8  J/s ‘m2* 
a  =  the  droplet  albedo  =  0.14 

p 

L  *  the  solar  insolation  rate  (J/m  *s) 
k  »  the  heat  transfer  coefficient  (J/m2,s‘K) 

■  kQ  (2  +  0.6  Re1//2Pr1//^)/d  (Reference  6) 

a 

where  k  =  the  thermal  conductivity  (J/m*s*K) 

Pr  =  the  Prandtl  number 
AH  =  the  latent  heat  of  vaporization  (J/kg) 

The  values  chosen  for  e&,  ed,  and  a  are  those  recommended  in 
Reference  B-8.  For  most  calculations  a  fairly  high  solar  rate  of 
1000  watts  per  square  meter  is  used  since  the  fuel  droplets  are 
released  at  high  altitude.  The  thermal  conductivity  and  Prandtl 
number  at  the  evaporating  film  are  essentially  those  of  air,  and 
these  were  approximated  from  tabular  data  by  the  expressions: 

k  =  0.024  (1  +  0.0034  (Tfi  -  273-15)) 

cl  a 

Pr  =  0.713  (1  -  4.6  x  10’4  (T  -  273-15)) 

cl 

From  inspection  of  the  heats  of  vaporization  for  a  range  of  indi¬ 
vidual  hydrocarbons  as  well  as  JP-4  itself,  it  was  determined 
that  in  the  temperature  range  of  interest  (-40°  to  +40°C),  aFv  can 
be  approximated  quite  well  by  the  linear  formula: 
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273.15)) 


AHy  -  3.  7  X  105  (  1  -  0.0013  (Td 

Once  Q  has  been  calculated,  the  amount,  T(K),  that  the  droplet 
cools  during  the  interval  is  simply: 

T  -  Q/md  Cp 

where  md  ■  the  droplet  mass  (kg) 

Cp  ■  the  droplet  heat  capacity  (J/kg/K) 

The  heat  capacity  of  petroleum  fuels  can  be  estimated  by 
(Reference  B-7) : 

Cp  -  (181  +  0.8  Td)/pd 

The  new  droplet  temperature,  mass,  composition,  and  height  are 
then  used  as  the  initial  conditions  for  the  next  interval.  This 
process  continues  until  the  droplet  has  accommedated  itself  to 
the  ambient  air  temperature.  The  steady-state  difference,  AT’ 
(K),  between  the  droplet  temperature  and  the  air  temperature  is 
determined  by  setting  Q  equal  to  zero  and  solving  for  AT'  using 
the  Newton-Raphson  iterative  method.  The  droplet  temperature 
thus  obtained  is  compared  with  the  droplet  temperature  which  was 
assumed  for  the  evaporation  calculations  and  if  the  discrepancy 
is  greater  than  0.5K  the  interval  is  rerun  using  the  new  value  of 
Td.  This  step  is  repeated  until  a  self-consistent  droplet  tem¬ 
perature  is  obtained. 

A  flow  chart  of  the  model  is  shown  in  Figure  B-2. 


APPENDIX  C 


FUEL  DROPLET  PROGRAM 

The  FORTRAN  IV  computer  program  developed  to  simulate  the  evapora¬ 
tion  and  free-fall  of  a  fuel  droplet  is  listed  in  this  appendix. 
The  variable  names  used  in  the  program  are  described  in  Table  C-l 
whioh  follows  the  listing. 
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PROGRAM  FRF^&l.  (  I  NPUT *  OUT  PUT  .  T  APF  7  =  I NPUT  ,  T  APFMOUTPUT  > 

SIMULATES  T  nF  FREE  FALL  AnU  EVAPORATION  OF  FUEL  DROPLETS 
REAL  m<40)  «Mn,MKv(40)  *M[  ,m,h.f.S*MOLF  (40)  *MR<40)  »MREL  *MS«MSAVE  (40)  * 
gM>TF P»MWI4«)  .NUM,AilJH,NAMl  ,NAM?,NAM3 
dimension  OFMS(40)  t«cv  (40)  , TR<40) ,VC<40> ,K(R) ,PCT<50)  tPCR(40i  * 

2-'  -Ni  (9) 

SET  CONSTANTS 
NMAX  = 10000 
TMAX=600000. 

RMMINs.OOOS 
OMMINe, 005 
DM*AX=.Q5 
OHMlNalO. 

DHMAX=1 00« 

EA=.75 
ED»,95 
ALB  =  . 14 
RS=1000. 

INPUT  COMPOSITION 
I  READ ( 7* 1 0 )  NCOMP 
IF (EOF (7) )  330*105 
i  WRITE (6*5) 

DO  110  1=1* NCOMP 

RE AO ( 7  *  20 )  NAM1,NAM2«NAM3.PCV(I) iMK(I)*TB(I) ♦ DENS (I) 

WRITE (6,30)  I,NAM1,NAM2,NAM3,PCV(I) fMW(I) * TB ( I > ♦ DENS ( I ) 

TR(I)=TR(I>+?73.15 

DENS ( I ) = 1 000 • *DENS ( I ) 

OCFB=l./ ( 1.  +  . 001* ( TB  < I » -293.15) > 

I  VC<I)=MWm/(DENS(I)*DCFB) 

READ ( 7  *  25 ) K 

K  CONTAINS  THE  LIST  OF  COMPONENT  NUMBERS  TO  BE  PRINTED  OUT 
READ (7*31) INn 

I ND=  1  FOR  SINGLE  DROP.  0  FOR  DISTRIBUTION 
IF(INO.EO.I)  GO  TO  122 

INPUT  DISTRIBUTION  (MULTIPLE  DROPLET  CASE) 
i  REAn(7,32)NR.DD*OX 
IF ( EOF ( 7 ) )  330  *  1 1 6 
.  REAO(7,34) (PCT(I) ,I=1*NB) 

DU=DD 

WRITE (6*36) 

DO  llfl  1=1. NR 
WRITF (6»38)DD»PCT (I) 

I  D0=00+DX 
DMAX=DD-OX 
DMIC=OU 
DS A VE=DU 

INPUT  CASE  TO  BE  PUN  (MULTIPLE  OROPLET) 

I  READ(7,40)HI,HG,TG,V 
IF (EOF ( 7) )  330,121 
IFU'.EQ.O.)  V=175. 

WRITE(6.50)DX,HI,HQ,T6,V 
WRITE  <6*51 ) 

WRITE (6*55) K 
WRITE (6,56) 

GO  TO  130 

INPUT  CASE  TO  HE  RUN  (SINGLE  DROPLET) 

!  READ ( 7 , 42 )  DMIC.HI ,HG, TG, V 
IF (EOF ( 7) )  330,123 
I  WRI TE ( 6  *  44 )  DMIC.HI ,HG* TG 
WRITE (6*46)  K 
WRITE (6.48) 

TMAX=60000. 

initialize  variables 

I  T0=TG*.0065*HG»2I3. 15 
IF(V.EO.O.) V=175. 

J=  1 

TOT=0. 

DO  132  1=1, NCOMP 
»  PCR ( I ) =0 . 

’  DI=.000001*DMIC 
D=D  I 
H-h  j 

T=HSTEP=MSTEP=MI=UT=0* 

TSAVE  =  - 1  * 

IDONE=0 


TSTEP* 1 • 

NPAS»i 

TA*T0-.00h5*H 
CSs20,*Sim  <TA> 

DTSAV*-TA*V*V/CS/CS/5. 

DT*DT5AV 

ins*o 

BEGINNING  OF  CALCULATION  LOOP 

TA*T0-.0065*H 

PA=101315.* ( TA/TO) **5.256 

DENSA*.0034R4*PA/TA 

VISCA=1.458E-6/< 110.4*TA)*TA**1.5 

TD*TA-DT 

TF*.5*(TA*TD) 

DCFD=1./(1.*.001*<TD-293.15>  > 

ON  INTIAL  PASS  CALCULATE  DROPLET  VOLUME  AND  MASS 
IF (NPAS.GT . 1 )  GO  TO  160 
VO*. 5236*0**3 
MIxO. 

DO  150  1*1, NCOMP 

M ( I ) *DENS ( I ) *DCFO*PCV (I)*VD/100. 

MSAVE ( I ) *M ( t ) 

M I =MI ,M  ( I ) 

MD*MI 

DENS* *MO/VD 

CALCULATE  COMPOSITION  AND  SIZE  OF  DROPLET 
MOLESxO. 

DO  170  1*1 .NCOMP 
MOLES*MOLES*m  < I ) /MW  < I ) 

VD=0. 

DO  180  1*1, NCOMP 
VD=VD,M < I ) / (DENS ( I ) *OCFO> 

D=( 1.90986* VO)** ( 1./3. ) 

DMIC=1000000.*D 
IF(IOONE.GE.l)  GO  TO  186 
DO  185  1*1 .NCOMP 
MOLF(I)xM(I) /(MW(I)»M0LES) 

DENSD=MO/VD 

MREL=MD/MI 

00  190  1*1, NCOMP 

mr(I)=m{i)/msave<i> 

DETERMINE  WHETHER  TO  PRINT 
IF(INO.EQ.O)  GO  TO  200 
IF  < IDONE • GE . 1 )  GO  TO  218 
IF (MREL.GE.RMMIN)  GO  TO  197 
IDONE=3 
GO  TO  218 

IF(T.LE.TSAVE)  GO  TO  220 

TSAVE=T 

GO  TO  218 

IF ( IOONE.GF. 1 )  GO  TO  210 
IF(MREL.GE.RMMIN)  GO  TO  220 
IDONF.  =  3 

MULTIPLE  DROPLET  PRINTOUT 
WRITE (6,60) DU » DM IC,T/60. ,H,MREL 
DO  211  1*1,9 
KK  =  K  < I ) 

PMR( I)sMR(KK) 

WRI TE ( 6 , 6 1 ) ( PMR (I), 1*1, 9) 

IF ( IDONE. IT ,?)  GO  TO  213 
IFUMAX.LT. 600000.)  GO  TO  213 
MREL=0. 

DO  212  1*1, NCOMP 
MR ( I ) =0 • 

WRI TE< 6,65) PCT(J)*MREL 
IF ( I  DONE i EQ«  ?  >  GO  TO  325 

INCREMENT  DROPLET  SIZE  AND  UPDATE  TOTALS  (MULTIPL 

DIJ=DU*DX 

DMIC*DU 

00  215  1*1, NCOMP 

PCR  <  I ) *PCR  < I ) +MR ( I ) *PCT ( J) 

TOT*TOT*MREL*PCT  <J) 

J=J*  1 

IF  NOT  DONE,  RE- INIT I AL I ZE  AND  BEGIN  NEXT  DROPLET 
IF(DU.LE.DMAX)GO  TO  137 


(MULTIPLE  DROPLET  CASE) 
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IP  DONE*  PR l NT  TOTALS  AND  CHECK  FOR  ANOTHER  CASE 
WRITE (6,70) TOT 

216  WRITE (6, 75) (1,PCR(1)*PCV(I)*QENS(I) *DCFD/DENSI/100. ,I«1 tNCOMP) 
OU=DSAVE 
DMlCaOU 
GO  TO  120 

SINGLE  DROPLET  PRINTOUT 

218  WRITE (6, 80) T/60. «H,DMIC*UT»MREL 
00  219  1=1,9 

KK=K(I) 

219  PMR ( I ) =  MR  <KK ) 

WRITE (6,81) (PMR ( I ) ,1=1,9) 

WRITE (6,82) OT 

IF  DONE,  CHECK  FOR  ANOTHER  CASE  (SINGLE  DROPLET) 

IF(IOONE.GE.l)  GO  TO  325 
CALCULATE  FREE  FALL  AND  EVAPORATION 

220  Q=SQRT< (4.*nENSA*0ENSD*9.8*D**3)/3.)/VISCA 
QLN=AL0G(0) 

REY=EXP (-3. 13+2.06*GLN-.083*QLN*QLN> 

UT=VISCA*RFY/ (D*0ENSA> 

HSTEP=UT*TSTEP 

MSTEP=0. 

00  225  1=1 , NCOMP 

DIFY=2.66E-5*SQRT<1./MW(I)+.0345>*TF*TF/(PA*<VC(I>**<1./3.>*.31)** 

12) 

NUM=2.+.6*SQRT (REY) *<VISCA/(DENSA*DIFY> )**(l./3.) 

HH=NUM*DIFY*MW ( I ) /(8314,34*0*TA) 

PV=EXP (20.53-2899./ (385. 15*TD/TB(I >-62.3) ) 

MEV  ( I ) =3. 14 1 6*0*D*HM*PV*H0LF ( I ) *TSTEP 
MS=M£ V ( I ) 

IF (MS.GT.M ( I ) )  MS=M(I) 

225  MSTEP=MSTEP*“S 
NpAS=NPAS+ 1 

IF(NPAS.LE.NMAX)  GO  TO  230 
IOONE=2 

SKIP  HEAT  BALANCE  IF  CALCULATING  DROPLET  COOLING 

230  IF ( lOS.EQ. 1 ) GO  TO  234 

CALCULATE  STEADY-STATE  TEMPERATURE  DIFFERENCE 
OLDTsOT 

TKA=. 024+8. IF -5* (TA-273. 15) 
pR=.713-3.3E-4* (TA-273. 15) 

NUH=2.+.6*S0RT(REY)*PR**(l./3.) 

HH=TKA«NUH/0 

DH V AP=3.7E5* ( 1.-.0013*(TO-273.15>  > 

231  DTNW=5.67E-8*EA*T  A**4-5 • 67E-8*ED*TD**4+ .25* ( 1 .-ALB) *RS*HH*DT-DHVAP 
l*MSTEp/(3.1416*0«0*TSTF.P) 

TD=TD+0TNW/(HH+2.27E-7*E0*TD**3> 

IF  <A8S(TA-DT-TD) .LE..1)  GO  TO  232 
OT  =  T  A-TO 
GO  TO  231 

232  IF (ABS(TA-TO-OLDT) ,LE..5)GO  TO  234 
OT=(TA-TD+OLOT)/2. 

GO  TO  140 
234  DT  =  T A-TD 

SKIP  DROPLET  COOLING  AFTER  STEADY-STATE  TEMPERATURE  IS  REACHED 
IF  < IDS.E0.2) RO  TO  240 
IF ( IDS.EO. 1 ) GO  TO  238 
IF (OT.GT.OTSAV)GO  TO  236 

-  SET  INDICATOR  THAT  DROPLET  COOLING  IS  NO  LONGER  NEEDED 
IOS=2 
GO  TO  240 

RESTORE  DT  AND  SET  INDICATOR  TO  CALCULATE  DROPLET  COOLING 
236  IDS=1 

DT =DTSA V 
GO  TO  140 

CALCULATE  DROPLET  COOLING  AND  RESET  INDICATOR  TO  HEAT  BALANCE 
238  I DS=0 

Cp=4184.*< 181 „+.8*TD) «VO 
DHVAP=3.7E5« ( 1 . -. 00 1 3* < TO-27 3. 1 5 ) > 

OCOOL=(EA*5.67E-8*TA**4-ED*5.67E-B*TD**4+.25*U-ALB)*RS+TKA*NUH* 
20TS A V/D) *3. 1416*0*0*1 STEP 
DCOOL= (DCOOL-MSTEP*DHVAP) /CP 
DTSAV=DTSAV-OCOOL 
DCSA V=DCOOL 
DT=DTSAV 
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HALVE  INCREMFNT  UNTIL  WITHIN  LIMITS 
)  IF(HSTEP.GT.nHMAX)  GO  TO  250 
IF(MSTtP/HI,l.E.OMMAX)  GO  TO  275 

1  aiFsams;?: 

MSTEP*0. 

DO  260  1*1  , NCOMP 
MEVm=MEV(I)/2. 

MS*MEV(I) 

IF(MS.GT.M(I)>  MSsMd) 

)  MSTEP*MSTEP»MS 

IF ( I DS. 60, 2) GO  TO  240 
OCSAV*DCSAV/P. 

DTSAV*DTSAV»OCSAV 

OT=DTSAV 

UPDATE2VARIARLES  WITH  RESULTS  FOR  INCREMENT 
i  HS=H 

H=H-HSTEP 

IF(H.LT.HG)  GO  TO  310 
T=T*TSTEP 

IFCT.GE.TMAXJGO  TO  321 

MO=MO-MSTEP 

DO  280  1  =  1  *NCOMP 

IF(MEV(I) .GT.M(I) )  MEV ( I ) *M { I ) 

)  M ( I ) =M ( I ) -MEV ( I > 

INCREASE  DURATION  OF  NEXT  INCREMENT 
TSTEP*2.*TSTEP 

IF (mo/MI .GE.RHMIN)  GO  TO  285 

I00NE=3 

GO  TO  140 

i  IF(HSTEP.LT.nHMIN)  GO  TO  290 
IF(MSTEP/MI.LT.OMMIN)  GO  TO  290 
GO  TO  140 
)  TSTEP=2.*TSTFP 
GO  TO  140 

END  OF  CALCULATION  LOOP 
INTERPOLATE  °0INT  OF  GROUND  IMPACT 
)  FAC=(HS-HG)/HSTEP 
IOONE=l 

i  TSTEP=TSTEP*FAC 
HSTEP=HSTEP*FAC 
MSTEP*0. 

DO  320  1=1 ,NCOMP 

MEVU)=MEV(I)*FAC 

IF (MEV ( I >  «GT.M ( I ) )  MEV<I)=M<I) 

MSTEP=MSTEP*MEV(I) 

)  M(I)=M(I)-MEV(I) 

T=T*TSTEP 
H=HS-HSTEP 
MD=MO~MSTEP 
GO  TO  140 

INTERPOLATE  POINT  OF  TERMINATION  BY  TIME 
l  FAC=(TMAX+TSTEP-T)/TSTEP 
IDONE=3 
T=T-TSTEP 
GO  TO  315 
>  WRITE (6,85)  NP AS 
. IF(IND.EO.l)  GO  TO  122 
GO  TO  214 
)  WRITE (6*90) 

STOP 

j  FORMAT (*1*»//»47X,*  COMPOSITION  OF  SIMULATED  FUEL  MIXTURE  *»///»19 
PX»*  COMPONENT*, 35X ,* VOLUME* , 1  OX, *MOL* ♦ 1 0X,*8OILING*,7X, ‘DENSITY**/ 
3, POX,*  NUMBER*, 8X, *SPEC I ES*,22X,*PERCENT*,8X,*WEIGHT*»8X,*PT,  (Cl* 
4 , 8X , *G/CC* , / ) 

}  FORMAT i 1 2) 

)  FORMAT (2A10,A9,F5.2,5X,2(F5.1,5X) ,F5.2) 
j  FORMAT (912) 

)  FORMAT  < *  *,22X,I2,10X,2A10,A9,1X,F5.2,10X»2(F5.1»9X) ,F5.2) 
l  FORMAT  (ID 

»  FORMAT (I?»8X,F5,0,5X»F5.0) 

V  FORMAT (5 (F5> 1 ,5X) ) 

i  FORMAT (////, 86 X,*  DROPLET  DISTRIBUTION  *,///,56X»*  DI AMETER*»5X»*W 
2EIGHT*»/,56X«*  (MICRONS) *,4X,*PERCENT*»/> 


GO  TO  290 


GROUND  IMPACT 


MEV ( I ) =M ( I ) 


38  FORMAT**  **Sox*F«>.0,6X*FS.  1) 

'%H8SSr.<:i;5^?.S^V55.5IgSS^vAE??TAS«J*IifiT-lMJi^oirKeB  m 

3X,*  LOCAL  GROUND  LEVEL  (METFRS)  !*.  3X.F#,  n^l?SIERS)  8  *3x*”6‘°*/l® 
*6  FORMAT (***«6?X*9(I2»4X) ) 

f^s.v^cT?svK>srft.:rsK5*iiMi&f*«i&J||J?®sar«i!ig:ft 

56  FORMAT(****il6X,*INITlAL  MASS**/) 

w-sgcs?  •«« 

FORMST(*.i:sflX^9(2X^F4j2n’r5*0’7X*F6*4*6X»F6-4»4X*9‘2X»f6.2)> 

p Ariii  *  *■  /  j*  -  >  . .  ..  __  7  •  L  '  * 


80 

81 

M2 

85 

90 


FORMAT (***«1?3X*F5.1) 

FORMAT (//,9X«*  TOTAL  NUMBER 
FORMAT**!  END  OF  JOB*) 


OF  ITERATIONS  «  **I7) 


1.02 


TABLE  0-1.  IDENTIFICATION  OF  PROGRAM  VARIBLES 


ALB  droplet  albedo 

CP  droplet  heat  capacity 

CS  local  speed  of  sound 

D  droplet  diameter 

DCFB  density  correction  factor  at  boiling  point 

DCFD  density  correction  factor  at  droplet  temperature 

DCOOL  droplet  cooling  during  interval 

DCSAVE  same  as  DCOOL,  for  use  during  interval  halving 

DD  droplet  diameter  in  microns  (during  input  of  distribution) 

DENS(I)  density  of  component  I 

DENSA  density  of  air 

DENSD  droplet  density 

DENSI  initial  droplet  density 

DHMAX  maximum  droplet  fall  during  interval 

DHMIN  minimum  droplet  fall  during  interval 

DHVAP  heat  of  vaporization 

DI  initial  droplet  diameter 

DIFY  diffusivity 

DMAX  diameter  of  largest  droplet  in  distribution 

DMIC  droplet  diameter  in  microns 

DMMAX  maximum  evaporation  during  interval 

DMMIN  minimum  evaporataion  during  interval 


DSAVE  diameter  of  smallest  droplet  in  distribution 
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DT 


difference  between  air  and  droplet  temperature 
DTNW  new  guess  for  steady  state  droplet  temperatuare 
DTSAV  saved  value  of  DT 

DU  initial  droplet  diameter  in  microns 

DX  diameter  Interval  for  distribution  of  droplets 

EA  emissivity  of  air 

ED  emissivity  of  droplet 

PAC  fraction  of  interval  before  droplet  reached  the  ground 

H  droplet  altitude 

HG  altitude  at  ground 

HH  heat  transfer  coefficient 

HI  initial  droplet  altitude 

HM  mass  transfer  coefficient 

HS  saved  value  of  H 

HSTEP  change  in  altitude  during  interval 

I DONE  flag  for  completion  of  calculation 

IDS  flag  for  droplet  cooling  calculation 

IND  flag  for  single  droplet  versus  distribution 

K(I)  number  of  component  to  be  printed  out  in  I'th  position 

KK  same  as  K(I)  (during  printout) 

M(I)  mass  of  component  I 

MD  mass  of  droplet 

MEV(I)  mass  evaporated  during  interval  (of  component  I) 

MI  initial  mass  of  droplet 

MOLES  total  number  of  moles  in  droplet 

MOLP(I)  mole  fraction  of  component  I 


MR(I)  fraction  of  mass  remaining  of  component  I 

MREL  fraction  of  mass  remaining  of  droplet 

MS  saved  value  of  MI> 

MSAVE(I)  saved  value  of  M(I) 

MSTEP  mass  evaporated  during  Interval  (of  droplet) 

MW(I)  molecular  weight  of  component  I 

NAM1  name  of  component  (during  input  and  printout) 

NAM2  name  of  component  (during  input  and  printout) 

NAM3  name  of  component  (during  input  and  printout) 

NB  number  of  bins  in  droplet  distribution 

NCOMP  number  of  components 

NMAX  maximum  number  of  Iterations 

NPAS  number  of  iterations 

NUH  Nusselt  number  for  heat  transfer 

NUM  Nusselt  number  for  mass  transfer 

OLDT  previous  value  of  DT 

PA  local  air  pressure 

PCR(I)  fraction  of  mass  remaining  in  droplet  distribution  bin  I 
PCT(J)  percent  mass  in  droplet  distribution  bin  J 
PCV(I)  percent  by  volume  of  component  I 

PMR(l)  fraction  of  mass  remaining  of  component  I  (during  printout) 
PR  Prandtl  number 

PV  vapor  pressure 

Q  variable  used  to  calculate  Reynolds  number 

QLN  log  of  Q 

REY  Reynolds  number 
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RMMIN  minimum  fractional  mass  remaining  of  droplet 
RS  solar  insolation  rate 

T  time 

TA  local  air  temperature 

TB(I)  normal  boiling  point  of  component  I 

TD  droplet  temperatuare 

TP  temperature  of  evaporating  film 

TG  temperature  at  ground-level 

TKA  thermal  conductivity  of  air 

TMAX  maximum  time  permitted 

TOT  total  mass  reaching  ground 

TSAVE  saved  value  of  time 

TSTEP  duration  of  Interval 

TO  temperature  at  zero  altitude 

UT  terminal  velocity  of  droplet 

V  true  airspeed  of  aircraft 

VC (I)  molar  volume  of  component  I  at  its  normal  boiling  point 

VD  droplet  volume 

VISCA  local  viscosity  of  air 


iOf> 


APPENDIX  D 


SAMPLE  FUEL  DROPLET  PROGRAM  RESULTS 

This  appendix  presents  the  predictions  of  the  fuel  droplet  evapor- 
tion  and  free-fall  program  for  five  cases  involving  a  single 
droplet  and  for  five  cases  involving  a  distribution  of  droplets. 
The  assumed  initial  droplet  composition  and,  for  the  latter 
cases,  the  droplet  size  distribution  is  shown  before  the  first 
case. 
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INITIAL  DISTRIBUTION 


Hq  TAC/SGPA 
Hq  SAC/SGPA 
Hq  USAPE/SGB 
Hq  PACAP/SGPE 
Hq  AAC/SGB 
Hq  APLC/SGB 
PAA/AEQ-10 
Hq  TAC/DEEV 
Hq  USAFE/DEPV 
AMRL/CC 
USAFSAM/CC 
ASD/CC 
AFOSR/CC 
AEDC/CC 

USAFRCE/WR/DEEV 

USAFRCE/CR/DEEV 

USAPRCE/ER/DEEV 

DTIC/DDA 

Hq  AFSC/SGB 

NAPC/Code  PE  71  AFK 

Hq  AFSC/dlwm 

Hq  AFSC/sdne 

Hq  USAF/LEEV 

0SAF/M7.Q 

OSAF/OI 

AFIT/L3GM 

AFIT/Library 

AFIT/DE 

R&D/EQ/Code  3021 
OEHI./CC 
Hq  AFESC/DEV 
USAFSAM/EDE 
Hq  AFISC 
AUL/LSE  71-249 
HqUSAFA/Library 
Hq  AFESC/TST 
OL-AD/OEHL 
OUSDR&E 
Hq  AAC/DEV 
Hq  AFLC/DEPV 
Hq  USAF/SGES 
EPA/ORD 
AMD/RDU 
Hq  AFSC/SGPA 
Hq  USAF/LEEVP 
Hq  USAF/LGYF 
AFIT/DEM 
Hq  ATC/DEEV 


(The  reverse 


1  U3AF  Hosp/SGB 
1  R&D/EQ/DARD-ARE-E 
1  Hq  AFESC/RDVC 
1  AMD/RDB 
1  Hq  AFSC/DEV 

1  USAFSAM/VNL 

2  AFGL/LKD 
1  ASD/DEP 
1  NEPSS 

1  Hq  SAC/DEV 
1  Hq  PACAF/DEEV 
1  AMRL/THE 
1  AFAPL/SFF 
1  AFOSR/N 
1  AEDC/DOTR 
1  Hq  MAC/DEEV 
1  AFOSR 

12  Hq  USAFE/DEVS 
1  Hq  MAC/SGPE 
1  23  CES/DEEV 
1  USCG  ( G-WEP-1/7 3 ) 

1  Hq  AFESC/RD 
1  Hq  AFESC/RDVA 
1  Army  Environmental 
1  Hygiene  Agency-HSE-EA 
1  OASD/( I&L)EES 
1  ARPA 

1  AFMSC/SGFA 

1  Hq  AFRES/DE 
^  EPA/ESRL 

4  O’ CGD9  (MEP) 

2  Hq  AFESC/RDV 
Hq  AFESC/WE 

j-  AFATL/DLODL 
1  AFWL/'SUL  (Tech  Lib) 

1  AFTEC/SGB 
1  Hq  AFRES/SGB 
1  ^TFW/DOV 
1  Hq  AFESC/RDVCA 
1  1  Med  Svc  Wg/SGB 
1  NAVFAC/Code  111 
1  Cnem  Abstracts  Ser 
1  NCEL/Code  15111 
1  USCG/GDD 
1  EPA/Corvallis 
1  EPA/Athens 
1  Hq  ATC/SGPAP 
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